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THE VARIABLE COEFFICIENT THIN OBSTACLE PROBLEM: 
CARLEMAN INEQUALITIES 

HERBERT KOCH, ANGKANA RULAND, AND WENHUI SHI 


Abstract. In this article we present a new strategy of addressing the (vari¬ 
able coefficient) thin obstacle problem. Our approach is based on a (variable 
coefficient) Carleman estimate. This yields semi-continuity of the vanishing 
order, lower and uniform upper growth bounds of solutions and sufficient com¬ 
pactness properties in order to carry out a blow-up procedure. Moreover, the 
Carleman estimate implies the existence of homogeneous blow-up limits along 
certain sequences and ultimately leads to an almost optimal regularity state¬ 
ment. As it is a very robust tool, it allows us to consider the problem in the 
setting of Sobolev metrics, i.e. the coefficients are only W 1,p regular for some 
p > n + 1. 

These results provide the basis for our further analysis of the free boundary, 
the optimal (C 1,1 / 2 -) regularity of solutions and a first order asymptotic ex¬ 
pansion of solutions at the regular free boundary which is carried out in a 
follow-up article, [KRS15], in the framework of VE 1,P , p > 2(n + 1), regular 
coefficients. 
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1. Introduction 

In this article we present a new, very robust strategy of analyzing solutions 
of the (variable coefficient) thin obstacle or Signorini problem. Let a 1 - 7 : Bf —> 
be a symmetric, uniformly elliptic tensor field which is W 1,p , p > n+ 1, 
regular and let B+ := {x £ B\ C R” +1 | x n +i > 0} denote the upper half-ball. Then 
consider local minimizers of the constrained Dirichlet energy: 

J(w) = j a 1 -* (diiv)(djw)dx, 

Bt 

where we use the Einstein summation convention and assume that 

w £ JC := {v £ v > 0 on B[ := B+ fl {x n +i = 0}}. 

Thus, local minimizers of this variational problem solve a uniformly elliptic diver¬ 
gence form equation in the interior of the upper half-ball, on this set they are “free”. 
However, on the codimension one surface B[ they obey the convex constraint w > 0. 
In this sense the obstacle is “thin”. 

Due to results of Caffarelli [Caf79], Kinderlehrer [Kin81] (who work in the setting 
with C 1 ’ 7 coefficients) and Uraltseva [Ura87], local minimizers are C 1,a regular for 
some a £ (0,1/2] and satisfy a second order elliptic equation with Signorini (or 
complementary ) boundary conditions: 

dia^djW = 0 in Bt , 

( 1 ) 

w > 0, —a n+ ’ 0 djW > 0 and w(a n+1 ’^djw) = 0 on B[. 

In the sequel we study (1) with the aim of obtaining optimal regularity estimates 
for its solutions as well as a better understanding of the free boundary, i.e. the set 
r„, = 3b'{ x £ B[ | w(x) > 0}, which separates the coincidence set, A„, := {x £ 
B[ | w(x) = 0}, in which the solution coincides with the obstacle w = 0, from the 
positivity set, il w := {x £ B[\ w(x) > 0}. 

1.1. Main results. In studying the variable coefficient thin obstacle problem (1) 
in a low regularity framework (we only assume that the coefficients a lJ are in an 
appropriate Sobolev class), we introduce Carleman estimates as a key ingredient of 
obtaining first information on the solutions. Here the Carleman estimate replaces 
a variable coefficient frequency function approach. 

Carleman estimates have the advantage of being very flexible with respect to 
perturbations: After deriving a constant coefficient Carleman inequality, it is of¬ 
ten possible to deduce a variable coefficient analogue by perturbative techniques. 
Comparing the Carleman estimate with monotonicity arguments, we note that 
monotonicity of the frequency function (for harmonic functions) is equivalent to 
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logarithmic convexity of the L 2 norm. It depends on rigid calculations. In con¬ 
trast, Carleman inequalities provide relaxed convexity statements (c.f. Corollary 
3.1). The conjugated equation is central (c.f. (18)), and one stays in the more 
flexible context of PDEs. 

Making use of the Carleman estimate, we obtain the upper semi-continuity of 
the vanishing order (c.f. Proposition 4.2), lower and uniform upper growth bounds 
for solutions (c.f. Corollary 4.1 and Lemma 4.1), compactness of L 2 normalized 
blow-up sequences (c.f. Proposition 4.3). Moreover, we obtain that along certain 
sequences, the blow-up limits are homogeneous (c.f. Proposition 4.5). This in 
particular allows us to classify the blow-up limits when the vanishing order is less 
than 2. Combining this information then allows us to prove the following almost 
optimal regularity result: 

Proposition 1.1 (Almost optimal regularity). Let a u : B+ — > R< i ym l ' lx( ' n+I ' 1 be a 
uniformly elliptic, symmetric W 1,p tensor field with p £ (n + l,oo]. Assume that 
w is a solution of the variable coefficient thin obstacle problem (1). Then, with 
7=1-^ 

|Vw(t) — Vw(y)\ 

I C(t) IMIl 2 (b+) \ x ~ vV forpG (n+l, 2 (n+l)), 

~\ C\\M\ L 2 {B +)\x-y\ 1,2 \n(\x-y\) 2 for p £ [2(n + 1), oo], 

for all x,y £ Bf. Apart from the dependence on 7 for the first case, the constants 
2 

C are also functions of ||a lJ Hu/i,p(s+) ’P’ n - 

In particular, this improves the regularity estimates of Caffarelli [Caf79], Kinder- 
lehrer [Kin81] and Uraltseva [Ura87], by coming logarithmically close to the ex¬ 
pected optimal threshold of C 1,1 / 2 regularity if p > 2 (n + 1) (or even attaining the 
optimal C 1,7 regularity if p £ (n + 1, 2 (n + 1))). 

The Carleman estimate however permits us to obtain more information: The 
existence of homogeneous blow-up solutions directly allows us to classify the lowest 
possible vanishing rate at free boundary points without invoking the Friedland- 
Hayman inequality [FH76] (c.f. the proof of Proposition 4.6). Since blow-up solu¬ 
tions satisfy the constant coefficient thin obstacle problem, it is not surprising that 
this lowest blow-up homogeneity coincides with the lowest possible homogeneity of 
solutions of the constant coefficient thin obstacle problem, k = 3/2. As in the case 
of the constant coefficient problem, we show that there is a gap to the next possible 
blow-up homogeneity which is, k = 2. 

In our follow-up article [KRS15], we use this and the upper semi-continuity of 
the mapping r ro 9 x 1-9 k x to separate the free boundary into 

= r 3 / 2 (w) U (J P«(w), 

K> 2 

where P 3 / 2 (r 0 is a relatively open set of the free boundary, the so-called regular free 
boundary and (J r K {w) consists of all free boundary points which have a higher 

k>2 

order of vanishing. Working in the framework of fT 1,p metrics with p £ (2(n+l), 00 ], 
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we prove that the regular free boundary is locally a C 1,a graph for some a € (0, 1) 
in [KRS15]. This then allows us to improve the almost optimal regularity result 
from Proposition 1.1 to an optimal regularity result. Moreover, we identify the 
leading term in the asymptotic expansion of solutions of ( 1 ) at the regular free 
boundary and provide explicit error bounds. 

1.2. Context and literature. Let us comment on the context of our problem: 
Apart from the (constant coefficient) two-dimensional problem, which was com¬ 
pletely solved by [Lew72] and [Ric78], and despite various impressive results on the 
higher dimensional problem [Caf79], [Fre77], [Kin81], [Ura87], the optimal regular¬ 
ity of the thin obstacle problem was only resolved relatively recently by Caffarelli et 
al. [AC06], [ACS08]. In these seminal papers a frequency function was introduced 
as the key tool in studying solutions of the thin obstacle problem. 

Following this there has been great progress in various directions for the constant 
coefficient problem: Relying on frequency function methods, the ground breaking 
papers of Caffarelli et al. [AC06], [ACS08] establish the optimal regularity of so¬ 
lutions, as well as the C 1,a regularity of the regular free boundary. This has been 
further extended to the related obstacle problem for the fractional Laplacian [Sil07], 
[CSS08] and parabolic analogues [DGPT13]. Relying neither on frequency functions 
nor on comparison arguments, Andersson [Andl3] has shown that similar results 
also hold for the full Lame system. Recently, Koch, Petrosyan and Shi [KPS14] as 
well as De Silva and Savin [DSS14] have proved smoothness ([KPS 14] proved an- 
alyticity) of the regular free boundary. Moreover, Garofalo and Petrosyan [GP09] 
give a structure theorem for the singular set of the thin obstacle problem. 

The variable coefficient thin obstacle problems is much less understood. Here 
the best regularity result in the literature in the setting of Sobolev regularity is 
given by Uraltseva’s C 1,a regularity result: In [Ura87] she proves that for the 
thin obstacle problem with a kF 1,p , p > n + 1 metric, a 1 - 7 , there exists a Holder 
coefficient a £ (0,1/2], depending only on the Sobolev exponent p and the ellipticity 
constants of a iJ , such that the corresponding solution is in C 1,a (c.f. also [Ura89]). 
Recently there has been important progress in deriving an improved understanding 
in the low regularity setting: In [Gui09] Guillen derives optimal regularity results 
for solutions of the variable coefficient thin obstacle problem. He works in the 
setting of G 1,7 , 7 > 0 coefficients. This was generalized by Garofalo and Smit 
Vega Garcia [GSVG14] to Lipschitz continuous coefficients by using the frequency 
function approach. In a recent work in progress Garofalo, Petrosyan and Smit 
Vega Garcia [GPSVG] further extend this result and prove Holder continuity of the 
regular free boundary. 

1.3. Difficulties and strategy. Let us elaborate a bit further on the main diffi¬ 
culties in investigating the variable coefficient thin obstacle problem: The central 
problem that has to be overcome and that is reflected in all stages of our arguments 
is the low regularity of the metric. This requires robust tools. 

In this first part of our discussion of the variable coefficient thin obstacle problem, 
Carleman estimates, which are well-known from unique continuation and the study 
of inverse problems [Car39], [JK85], [KT01], [Isa04], [Rull4b], are used to handle 
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this low regularity setting: Although they are usually employed in the setting of 
Lipschitz metrics (as for instance the unique continuation principle in general fails 
for less regular coefficients), it is possible to extend these to the low regularity 
framework we are interested in. 

This allows us to carry out a blow-up argument and exploit information from 
the constant coefficient setting. For this we argue in two steps: 

• A Carleman estimate for solutions of the variable coefficient thin obstacle 
problem (f). This yields sufficiently strong compactness properties in order 
to carry out a blow-up procedure at the free boundary points. In particular, 
doubling inequalities (c.f. Proposition 4.3) are immediate consequences of 
the Carleman estimate. 

• A blow-up procedure. The good compactness properties deduced in the 
previous step permit to carry out a blow-up procedure with a non-trivial 
limit satisfying a constant coefficient thin obstacle problem (c.f. Proposition 

4.4). Thus, (for an appropriately normalized metric a lJ ) the blown-up 
solution is of the form: 

Aw = 0 in Bf, 

( 2 ) . 

w > 0, — d n +iw > 0 and wd n +iw = 0 on B[. 

However, the Carleman estimate yields further information: It is possible 
to show that the there are blow-up sequences such that the blow-up limits 
are homogeneous solutions of (2). Moreover, the homogeneity is given by 
the order of vanishing (c.f. Proposition 4.5). For a solution of the variable 
coefficient thin obstacle problem this then allows to exploit the existing in¬ 
formation on solutions of the constant coefficient thin obstacle problem and 
thus, for instance, obtain an almost optimal regularity result (c.f. Proposi¬ 
tion 4.8). 

1.4. Organization of the paper. Let us finally comment on the structure of the 
remainder of the article: In the next section, we first briefly recall auxiliary results 
(c.f. Proposition 2.1) and then introduce our notational conventions in Section 2.2. 
Following this, Section 3 is dedicated to our Carleman estimate, Proposition 3.1. In 
particular, we derive an important corollary, Corollary 3.1, from it. In Section 4 we 
then deduce crucial consequences of the Carleman estimate: We derive compactness 
properties (c.f. Proposition 4.3), carry out a blow-up procedure (c.f. Proposition 

4.4) and prove the existence of homogeneous of blow-up solutions (c.f. Proposition 

4.5) . This is then applied in proving the first (almost optimal) regularity result, 
Proposition 4.8. 

In our second article dealing with the thin obstacle problem we then use these 
results to further analyze the regular free boundary, to derive optimal regularity 
estimates and to obtain a first order expansion of solutions at the regular free 
boundary. 


2. Preliminaries 
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2.1. Auxiliary results and assumptions. In the sequel we recall certain auxil¬ 
iary results which allow us to reformulate and simplify the problem ( 1 ) in a particu¬ 
larly useful way. Moreover, we collect all our assumptions on the involved quantities. 

We start by recalling a (slight modification of a) result due to Uraltseva [Ura89], 
p.1183 which allows to simplify the complementary boundary conditions: 

Proposition 2.1. Let a lJ : B+ —> Rsym 1 ' X ^" +1 ' ) be a uniformly elliptic W 1,p , 

p > n + 1, tensor field and w £ H l (B+) be a solution to (1). Then for each point 

x £ B'i there exist a neighborhood, U(x), and a W 2,p diffeomorphism 
2 

T : U{x) fl Bf Abn B +, x i-A T(x) =: y, 

such that in the new coordinates w{y) := w(T~ 1 y ) (weakly) solves a new thin 
obstacle problem with 

d k b M (y)d t w = 0 in V n B(~, 
dn+iw <0, w > 0, w(d„+iu)) = 0 on V fl B [, 
where B(y ) = ( b kt (y)) = \det(DT(x))\- 1 (DT(x)) t A(x)DT(x)\ x=T _ 1{y) satisfies 

b n+1,t _ q Qn j or a n £ g |i j ... ) n}, 

and remains a uniformly elliptic W 1,p tensor field. 

Here the boundary data in (1) and (3) are interpreted in a distributional sense: 
It is known that one can see from the equation by local arguments that, if w £ H 1 is 
a weak solution in the interior, then a n+1 ’^djW | Xn+1 =o £ H -1 / 2 . This remains true 
in the new coordinates: b n+1 ’ n+1 d n+ \w\ Xri+l ^o £ If -1 / 2 . Since b n+1 ’ n+1 is bounded 
from below and in W 1,p , the multiplication by its inverse defines a bounded op¬ 
erator on H 1 ! 2 and l? -1 / 2 . Thus, wd n+ iw\ Xrl+1 =o defines a distribution on the 
boundary which we require to be zero. Vice versa: If w £ H 1 , dkb kl diiv = 0 in 
V fl B+ and d n+ \w < 0, w > 0 and wd n+ \w = 0, then w is a local solution to 
the thin obstacle problem. In any of the two formulations the regularity theory of 
Uraltseva [Ura87] applies. 

By a further (affine) change of coordinates we may, without loss of generality, 
assume that 

\b u (x)-5 u \ < C'l^P, 

for 7 = 1 — an d | 1 ^ >r ^ denotes the Kronecker delta. In the 

' p | 0 else, 

sequel, for convenience, we will always assume that we are in a sufficiently small 
coordinate patch such that our thin obstacle problem is formulated in this way. 
This allows us to exploit the boundary conditions very efficiently. So, without loss 
of generality, we pass from ( 1 ) to considering 


( 4 ) 


dia^djW = 0 in B +, 
w > 0 — <9„+i w > 0 and w(d n +iw) = 0 on B[, 


with 

(AO) |MIl2 ( b+(o)) = 
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(Al) a l ’ n+1 (x', 0) = 0 on R" x {0} for i = 1,..., n, 

(A2) a 13 is symmetric and uniformly elliptic with eigenvalues in the interval 
[ 1 / 2 , 2 ], 

(A3) a IJ G W 1,p (Bi (0)) for some p £ (n + 1, 00 ], 

(A4) a ij {0)=5 ij . 

Moreover, throughout the paper we assume that the solution, w, of (4) satisfies 
(W) there exists a £ (0, 5 ] such that w £ C 1,a . 

Assumption (AO) is a normalization, which is no restriction of generality, since 
positive multiples of solutions are solutions. Property (Al) may be assumed because 
of Proposition 2.1. By Morrey’s inequality and assumption (A3) 

|a^(x)-^|<C'||Va^|| iI)(Bi+( 0 )) |xr, 

where 7 = 1 — and hence the condition (A2) on the eigenvalues holds after 
rescaling if necessary. The regularity assumption (W) does not pose any restrictions 
on the class of solutions, as Uraltseva [Ura87] proves that this property is true for 
any H 1 solution of (4). We however stress that most of the paper is independent 
of this regularity assumption in the sense that assumption (W) is only invoked 
in Proposition 4.6. Apart from this, all the other necessary regularity results are 
proved “by hand” in the paper. 

2.2. Notation. In this subsection we briefly introduce some notation. We set: 

• R ” +1 := {x £ R” +1 | x n+1 > 0}, RK +1 := {x £ R n+1 | x n+l < 0}. 

• Let Xq = (£q, 0) G R+ +1 . For an upper half-ball of radius r > 0 around 
Xq we write B+{x 0 ) := {x £ R" +1 | \x — Xq\ < r}; the projection onto the 
boundary of R" +1 is correspondingly denoted by B' r [x 0 ) := {x £ R"j \x — 
Xq\ < r}. If xo = (0,0) we also write 5+ and B' r . Analogous conventions 
are used for balls in the lower half sphere: B~(x 0 ). Moreover, we use the 
notation B"(x 0 ) = B' r {x 0 ) fl {x n = 0}. 

• Annuli around a point xo = (^o,0) in the upper half-space with radii 0 < 
r < R < 00 as well as their projections onto the boundary of R ” +1 are 
denoted by by A+ R (x 0 ) := B+(x 0 ) \ B+(x 0 ) and A' r R (x 0 ) ~ B' R (x 0 ) \ 
B' r (x 0 ) respectively. For annuli around xq = (0,0) we also omit the center 
point. Furthermore, we set A r ^(x 0 ) := A)/ B (a:o) U A~ R (x 0 ). 

• f } w := {xeR”x {0} | w{x) > 0} denotes the positivity set. 

• r™ := Ob'^w is the free boundary. 

• A^, := B[\ Q w is the coincidence set which we also denote by A^. 

• We use the symbol A < B to denote that there exists an only dimension 
dependent constant, C = C(n ), such that A < C(n)B. Similar conventions 
are used for the symbol >. 

3. A Carleman Estimate 

3.1. A Carleman estimate and variations of it. In this section we introduce a 
central tool of our argument: The Carleman estimate from Proposition 3.1 allows 
us to obtain compactness properties for blow-up solutions (c.f. Proposition 4.4), to 
deduce the existence of homogeneous blow-up solutions (c.f. Proposition 4.5) and 
to derive the openness of the regular part of the free boundary (c.f. Proposition 
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4.2). It is a very flexible and robust tool replacing a variable coefficient frequency 
function argument. 


Proposition 3.1 (Variable coefficient Carleman estimate). Let n > 2 and 0 < p < 
r < 1. Let a 13 : A^ r R”J^ be a tensor field which satisfies 

(i) of 3 £ W 1 ' n+1 {A+ r ), 

(ii) the off-diagonal assumption (Al), 

(iii) the uniform ellipticity assumption (A2), 

(iv) the following smallness condition: There exists S = S(n) > 0 such that 

(5) sup || Va^|| A + < <5. 

p<r<r r,Zr' 


Assume that w £ H 1 (B+) with supp(w) C A^ r is a weak solution of the divergence 
form equation 

dia l3 djW = f in A+ r , 
w > 0, d n +\w < 0, wd n +iw = 0 on A' r , 


( 6 ) 


where f : A+ r —> R is an in A^ r compactly supported L 2 (A+ r ) function. Let <fi be 
the following radial weight function: 


4>(x) = </>(ln |x|) with (j) : 


<p(t) = —t + Co [ t arctan t — — ln(l + t 2 ) 


where 0 < Co 1 is an arbitrarily small but fixed constant. Then for any 7 £ (0,1) 
and any r > 1 we have 

( 7 ) 


e T ^|x| _1 (l + ln(|x|) 2 ) _ A 


+ r 2 


e T *(l + \n{\x\) 2 )~?X7w 


< c 0 

where 

( 8 ) 


L 2 (A+ r ) 

l C{n) (r 2 C(a y ) ||e T 0 |a:| 7 “ 1 w|| i2(j4 + r) + ||e r< ^|x|/|| i2(A + r) ) , 


L 2 (A+ r ) 


C(a 13 ) = sup 


A + 
A PA 


T 3 {x ) - S i3 


\xV 


+ sup |||x| 7 Vo* J '| 

p<r<r/2 


L"+i(A+ 2 -) • 


Remark 1 (Restrictions on r). In order to derive the Carleman estimate, it would 
not have been necessary to assume that r > 1. In fact the estimate is valid for 
arbitrary r > 0 if we replace the right hand side of (7) by 


T 2 

< 


e T ^|x| *(1 + ln(|x|) 2 ) 2 w 


- T 2 


+ ln(|x|) 2 ) 2 Vic 


L 2 (A+ r ) L 2 (A+ r ) 

c Q 1 C(n) (max{l, T 2 }C(a lJ ) ||e r 0 |x| 7 ~ 1 w|| i2(A + r) + ||e r 0 |x|/|| i2(A + r) 

+r^ e T< ^|x| _1 (l + lndd) 4 ) - ^ 

L 2 (A+ r ) 


Choosing 0 < r < 1 depending on r sufficiently small, then allows to absorb the 
error terms in this case as well (c.f. Remark 2). 


Remark 2 (Applications). In the sequel we will apply the Carleman estimate (7) 
with a metric a 13 such that a 13 (0) = 6 13 and a 13 £ W 1 ’ P (B^), p £ (n+ 1, 00 ]. Thus, 
by Morrey’s and Holder’s inequalities, for 7 = 1 — ^±1, C(a 13 ) is bounded by a 
constant which depends only on n and ||Va y ||^p(s+)- 
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In our applications, we will always consider r < Rq = Ro(r,n,p, || Va ‘ 3 Hip(s+)) 
sufficiently small, such that the smallness condition (5) is satisfied and moreover, 
the first term on the right hand side of (7) can be absorbed by the left hand side. 
More precisely, we consider i?o such that 

C{a ij )Rl < 5, t^Cq 1 C'(n)C'(a ij )i?g| Ini? 0 | < 1/4. 


Then the Carleman inequality (7) can be rewritten as 
(9) 


r 2 


e T *\x\-\l+\n(\x\) 2 )-^ 


L 2 (A+ r ) 


■ r 2 


e T?i (l + ln(|x| ) 2 )-5Vw 


L 2 (A+ r ) 


< c 0 1 C(n ) ||e T,i |x|/|| i2 ^ A+ j for 0 < p < r < R 0 and r > 1. 


Remark 3 (Carleman estimate for n = 1). In one dimension it is possible to prove 
a similar Carleman estimate as above. However, slight modifications are necessary. 
Instead of working with the W 1,2 semi-norm (which would be the scale invariant 
norm in (1 + l)-dimensions) of the metric and the smallness condition (5), we con¬ 
sider L p norms, ||Va® : '|| iP withp > 2. Moreover, all Sobolev embedding arguments 
are then replaced by interpolation inequalities of the form 

(10) |H| j. < c(p, n) || Vw ||£ 2 |M|* 2 p ■ 

L/ P~ z 


This then yields the following Carleman estimate 

( 11 ) 


3 

T 2 


< 


e T ^|x| *(1 + ln(|x|) 2 ) 2 w 


L 2 (A+ r 


- T 2 


+ ln(|x|) 2 ) 2 Vic 


L 2 (A+ r ) 


C 0 1 c(n,p) ( 'T 2 || Va lJ || iP ||e T ^|x | 7 1 w\\ L 2 {A + r) + || eT0 l a: l/||L2 (A + r . ) ) • 


Remark 4 (Generalizations). • With slight modifications, the Carleman es¬ 
timate is valid for more general (radial) pseudoconvex weight functions, 
(f>(x) = </>(ln(|x|)). In this more general setting (but with <fi(t ) ~ —t and 
\<p'(t)\ ~ 1 asymptotically as t —> — oo) the Carleman estimate reads 


T 2 


( 12 ) 


e r ^|x| 1 (j)' (cj)") 2 w 


L 2 (A+ r ) 


- T 2 




L 2 (A+ r ) 


< c 0 


l C(n) (r 2 C(a ij ) ||e r ^| 




I L 2 (A+ r ) 


5 I/| 


L 2 (A+ r ) 


)• 


where C(a ®1) is the constant from ( 8 ). 

• As we aim at absorbing the error terms on the right hand side of (7) or 
( 12 ) into the left hand side of the inequality if r > 1 is uniformly bounded 
by some To < oo, it is in general not possible to work with “linear” weight 
functions in ( 12 ). As a consequence, in the sequel we will always consider 
strictly pseudoconvex weight functions with the exception of Lemma 3.2. 

• In contrast to the setting of unique continuation, where the metrics under 
consideration are Lipschitz, our metric is only W 1,p , p > n + 1, regular. On 
the one hand, this requires additional care in the commutator estimates of 
the Carleman inequality. On the other hand, the low regularity of the metric 
also implies that we cannot hope for the unique continuation principle to 
hold in the setting of (4). Hence, the Carleman estimate can, for instance, 
not exclude free boundary points with infinite order of vanishing. On a 
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technical level these observations are manifested in the error contributions 
on the right hand side of the inequality (7) which carry higher powers of r 
than the analogous contributions on the left hand side. 

• It is possible to derive an analogous Carleman inequality for the setting of 
the interior thin obstacle problem and the case in which inhomogeneities 
are present in the thin obstacle problem. For a detailed discussion of this 
we refer to Section 5 in [ ARS15]. 

Before proving the Carleman estimate, we state H 2 estimates for the associated 
conjugated operator in dyadic annuli which will be proved in Section 3.3 at the end 
of this Section. 

Lemma 3.1. Let n > 2. Then there exists 5 = S(n) > 0 such that the following is 
true: Let a 13 : B+ —> be a uniformly elliptic, symmetric tensor field 

which satisfies (5). Further, assume that v : —> R with v £ H 1 (B(f) is a weak 

solution of 


dia^djV + t 2 {di4>)a^ {dj(f>)v — ra 13 (difydjV 

- T{dj4>)di{a l °v) - to , 13 {d^cjTjv = f in B+, 
d n +\v <0, v > 0, vd n +iv = 0 on B[, 
v = 0 on dBi, 

where 4>(x) = <^(ln(|ic|)) is defined as in Proposition 3.1 and f : B(~ —>■ R is an 
L 2 (B+) function. Then on each dyadic (half-)annulus A+ := B+_ m \B^_ m _ 1; 
m £ N, m > 1, and for each r > 1 

(14) ||V 2 u|| i2(A +) < T 2 || |a?| 2 HIa 2 (A+_ iU A+ U A+ + 1 ) + ll^lli 2 (A+_ 1 UA+UA+ +1 ) • 

Remark 5 (Modifications for n = 1). Analogous to the Carleman estimate, Lemma 
3.1 has a similar extension to n = 1. In that case, we replace the smallness condition 
(5) and all Sobolev embedding arguments by interpolation arguments as in (10). 

Proof of Proposition 3.1 . We interpret the equation as a perturbation of the Lapla- 
cian 

(15) Am = di ( a ij - 6 ij )d jW + f, 

and argue in two steps: First, we derive a Carleman inequality for the Laplacian 
with Signorini boundary conditions. In a second step, we explain how to treat the 
right hand side of (15). However, before coming to this, we recall that the boundary 
contributions in (6) are well-defined as distributions. 

Step 1: Regularity of w. In order to make sense of the boundary contributions, 
we refer to the discussion following (3). Moreover, using the H 2 regularity of w 
(c.f. Lemma 3.1), we may also deal with second derivatives of w in L 2 (B)() as well 
as L 2 gradient contributions on the boundary. An alternative method to deal with, 
for example, boundary contributions would have been a regularization mechanism 
(c.f. Remark 6) 
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Step 2: Carleman for the Laplacian with Signorini boundary conditions. We will 
show that 
(16) 


3 

r 2 

e T ^\x\ 1 (l + ln(x) 2 ) 2 w 

+ T 2 

e r ^(l + ln(|x|) 2 ) 2 ~Vw 



L 2 (A+ r ) 



< c 0 1 C(n) \\e T ^\x\Aw\\ L2{A + r) . 


As in the article of Koch and Tataru, [KT01], we carry out a change into radial 
conformal coordinates. 

To this end, we first consider the (n+ l)-dimensional Laplacian in polar coordinates 

77< 1 

d r + — d r + 

and then carry out a change into radial conformal coordinates. Setting r = e 4 , 
yields d r = e~ t dt , which in turn transforms the operator into 

e 2t (9 2 + (n — 1)9 1 + As«). 

Conjugating with e ~ t (which corresponds to setting w = e ~ t u) and mul¬ 
tiplying the operator with e 2t , results in an operator on R x S" which is of the 
form 

(17) £ + As „. 

In the sequel, we only work with this modified operator and prove a Carleman 
estimate for it. Under the described change of coordinates and using the comple¬ 
mentary boundary conditions combined with the H 2 regularity of w (which follows 
from Lemma 3.1), the original Cartesian boundary conditions 

w > 0, — d n+ \w > 0 and wd n+ iw = 0 on R", 


now read 


u > 0, v ■ Vgnu > 0 and u(y • V s n u) = 0 on R x S’ 

where v = (0,..., 0, —1) denotes the outward unit normal to 
In these new coordinates (16) can be phrased as 


n— 1 


'(1 +f 2 )“U 


< c 0 1 C(n) 


"b T 

L 2 ((—oo,0] x5p 
2 

+ A S n I U 


e^(l+t 2 )-W (t , S n)i 


L 2 (( —oo,0] xS?) 


,2 ( n - !) 


df- 


L 2 (( —oo,0] x5|) 

We conjugate the operator (17) with the only f-dependent weight, e r< ^ (and in 
particular, we define u = e T ^u). This leads to the following (up to boundary 
contributions) “symmetric” and “antisymmetric” parts of the operator: 


32 i _2 /q 7\2 ( n l ) 2 


(18) 


5 = 9 2 + r 2 (9 ^) 2 - 
A = -2 T(d t <j>)d t - rdlf). 


+ A S n, 


In order to derive the desired Carleman estimate, we expand the L 2 -norm of (A + 
S)u : 

\\(S + A)u\\ 

L 2 (M ) = II'S'mII i 2 (M) + \\ Au \\l 2 (m) +2{Su,Au) L 2( M) , 
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where M = R_ x 5”. In order to bound this from below, we have to control the 
last contribution 2(5u, Au)^ 2 ^ M y To this end, we integrate by parts and reformu¬ 
late it as a commutator. As our “symmetric” and “antisymmetric” operators are 
only symmetric and antisymmetric up to boundary contributions, this also leads to 
boundary integrals: 

(19) 2(Su, Au)l 2 (m) = {[S,A]u,u)l2 (m) + boundary contributions. 

A calculation shows that the boundary contributions are of the form 

( 20 ) 

— 8 t J d t <j){v -V s n u)d t ud'H Tl ^ 1 dt — 4r J d 2 (j){y -V s n ‘u)ud'H n ~ 1 dt, 

R _xS”-' R_xS '"- 1 

where v = (0,..., 0,-1) denotes the outer unit normal field associated with 5”. 
We begin by discussing these boundary contributions: Due to the radial dependence 
of (f> (which in Cartesian variables implies that all the terms involving d n +\<f> vanish 
on supp(ic) nB(), the integrands in (20) correspond to contributions of the form 


(x ■ Vw)d n +iw and wd n+ iU ), 


in Cartesian variables. By the complementary boundary conditions and the H 2 
regularity of w, both of these vanish on B[: In fact, w £ H l (A' p r ) and thus 
x ■ Vm = 0 almost everywhere where w{x) = 0 (c.f. for instance Theorem 6.17 
in [LL01]). In particular, x ■ Vic = 0 for almost every x € A w . Combined with 
the vanishing of d n +i w on fl w , this implies that the boundary contributions in the 
Carleman estimate vanish. 

We hence proceed to the bulk commutator contribution in (19): As the weight (j) is 
sufficiently pseudoconvex, the commutator is bounded from below. Indeed, 

([5, A]u,u) L 2 {M ) = 4T((d?(f))dtu,dtu) L 2 iM) + 4T 3 ((<9 t 2 ^)(<9 t ^) 2 u,ii) L 2 (M) 

- T{{d?4>)u,u) L 2 (M) . 


Using the explicit form of (j>, we have df<j)(t) < (9/8 )df<p(t) for all t, < 0. Hence, 
the last term can be absorbed into the other contributions if r > 1: 


([S’, A]u,u) L 2(m) > 4r ( dt(j)) 1,2 dtu 


L 2 (M) 


2 r d 


(d?4>) 1/2 (d t j))u 


L 2 (M) 


Hence, we obtain 


||(5 + A)u\\ l2(m) > \\S u \\l 2 (m) II^- u IIl 2 (m) 


+ 4r 


(d 2 4)) 1/2 d t u 


L 2 (M) 


2 T d 


(d?0) 1/2 (dt4>)u 


L 2 (M) 


Last but not least, we upgrade the gradient estimate from an estimate for the 
radial derivative to an estimate for the full gradient: Using the symmetric part of 
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the operator, we deduce 


(dt4>) 2 V S »M 


L 2 (M) 


{dUY 2 d t u 


L 2 (M ) 


= -((< 9 i 2 ^)u, A S n U ) L 2 (M) - {(d 2 4>)u,d 2 u) L 2 (M) 
(d 2 f))Y • V s »u)uM n - l dt 


/ 


IxS" 


( 21 ) 


= -({dt4>)u,Su) L 2 {M ) +T 2 {{dt^)(d t 4>) 2 u,u) L 2 (M) 
(n^l ) 2 


+ 


-({d t <t>)u,u) L 2(M) 


< r 


+ 


(d?(j>)u 
( n ~ l ) 5 


L 2 (M) 


-\\Su\\ 2 L 2 (m) + t 2 (d 2 4>)Hdt4>)u 


2 

L 2 (M) 


(d?j>Yu 


L 2 (M) 


Again we used the complementary boundary conditions to dispose of the boundary 
integral. Observing that all right hand side contributions in (21) can be absorbed 
into the left hand side of the Carleman inequality and that this remains true if 
we multiply (21) with the factor cr for a small, positive constant c, we obtain the 
desired full gradient estimate. 

Thus, inserting the changes we made in passing to our conformal change of vari¬ 
ables, i.e. w = e^~ t e~ T l > u, and recalling the changes in the volume element, results 
in the claimed inequality (16). 

Step 3: Bounds for the right hand side of (15). We now proceed to estimating 
the right hand side of (15) and (7). Applying the triangle inequality immediately 
leads to 

\\e T,l, \x\Aw\\ L2{A + r) < ||e r ^|a;|/|| i2(A +^ ) + \\e T(t, \x\ di(a ij - 6 ZJ )djW\\ L 2 {A + r) ■ 

Thus, it remains to control the second contribution. Setting v = e T l > w 1 in Cartesian 
coordinates after conjugation, we have 

e T *\x\d t {a ij - 5 ij )djW = \x\ [(dia^djV - + (a lj - 6 lj )d tj v 

- 2 r(a ij - + r 2 (a ij - - r(a ij - S'- 1 )(<),, o)i(. 


Since \a l l(x) — 5 lJ: | < sup 




\x \2 


|x| 7 < C(a l l) |x| 7 , we can directly bound all the 


contributions which neither contain derivatives of a lJ nor second order derivatives 
of v. More precisely, we have 


2T\\\x\{a ij 


L 2 


|®|(a y - 5 ij )(dij<t>)v I 


+ T 


L 2 


\x\(a^ - 


L 2 


< r ||(a u — 8 l °)djV || L 2 + r 

< TC(a l Y |||x| 7 3ju|| i2 +T 2 C(a lJ ) |||x 


|x| -'(a^ -8 ij )i 


I L 2 


■IT- 1 , 


lL 2 ’ 


where we used that |V<^(x)| < C |cc| —1 , \S7 2 (f(x)\ < C|a;| _2 . Thus, using that the 
gradient term can be controlled by the L 2 contribution on the right hand side of 
the last inequality (c.f. the comments at the beginning of the proof of Lemma 3.1), 
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these errors are of the form of the contributions on the right hand side of (7). 


We continue with the contribution which involves second derivatives of v. We 
first decompose the L 2 norm over B+ into L 2 norms on the dyadic (half) annuli 
A+:= B+_ m \B+_ m _ 1 -. 

IN ^ E 2_msu Pll%^ll^(A+) • 

meN At, 

On each annulus we invoke Lemma 3.1: 

ll^b' U llz, 2 (A+) ~ T 2 2 2m IMIi 2 (A+_ 1 UA+UA+ +1 ) + ll^lli 2 (A+_ 1 UA+UA+ +1 ) • 

Using the finite number of overlaps of the annuli in the right hand side contributions, 
we can sum over these. Combining this with the observation that 

sup \a ij (x) - 5 ij | < 2 and sup | a ij {x) - S ij \ < C(a ij ) 2“ m7 , 

At, A+ 

then results in 

\\\x\{a lJ - S lJ )di j v\\ L2{ B+) 

< 2 _m sup \a^{x) — | (r 2 2 2m \\v\\ L2(A+m _ iUA+mUA+m+i) 

m£ N A™. 

+ II/IIl 2 (A+_ 1 UA+UA+ +1 ) 


< 


Y (c{a ij )T 2 2 m M||«|| 

m<E N 
+ 2~ m 

< C(a l2 )T 2 |||a;| 7_1 u 


i 2 (-4+_iUA+UA+ +1 ) 

L 2 (A+_i u A+LIA+ + i)) 

+ ^ T |||®|/||j 


li 2 (B+)~ BI-^i j ili 2 (B+) ■ 

All of these are terms as in the right hand side of (7). 

Hence, it remains to consider the terms which carry derivatives on the coefficients 
<W . Introducing a partition of unity and using Holder’s inequality, we have 


\x\{dia^)djV 


l u \\ L 2 (B+) 


— 'y \ \\9iU 2 ||^ J n+ 1 (^+ ) || 2 ( n +H 


m£ N 


(B+) 


Here we assume that each of the functions r] m is supported in the annulus A 2 -m-i 
and satisfies derivative bounds of the form \djT] m \ < 2 m . Sobolev’s inequality then 
yields 


\x\r) m djV || 2 ( n+i) < 2 " 

L (B+) 


( 22 ) 


< 


\\VmdjV\\ . 2 1 "+| 1 ) ., + . 

|x|dijv|| 

L 2 (A+_iUA+UA+ +1 ) 
' H^' t ’lli 2 (A+_ 1 UA+UA+ +1 ) ' 


Invoking the H 2 estimates from Lemma 3.1 again, we obtain 

lll :r l%' i; llz, 2 ( J 4+) ~ T ^ IMIi 2 (A+_ 1 UA+UA+ +1 ) + ^ ll/lli 2 (A+_ 1 UA+UA+ +1 ) 

Using a combination of 

pi^L^HAt,) - s and ll 5iaZJ IL™+i(A+) - 2_mc V' j )> 
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we again sum these contributions over m: 


\x\{d 


i U \\L 2 (B+) 


s E ll». 

mG N 


& ^ r n+i(A + \ || |^|^7m^7^|| 2 ( n +i) 

UL> L ri-l (B+) 

~ \\di a ,? llz,'*+ 1 (A+) ( T 2 IMI_L 2 (A+_ 1 UA+UA+ +1 ) 


mGN 

2 _m I 


l_L 2 (A+_i u A+LJA+ + i )) 


< 


J2(C(a ij )r 2 2 r 




m(EN 


Ia 2 (a+_ 1 ua+ua+ +1 ) 


+ 52 m II/IIl 2 (a+_ iU a+ u a+ +1 )) 

< C(a lJ )r 2 1|l a: r 1+7 «|| I , 2(B + ) +<5 |IM/IIl+b+ ) 


This can be absorbed into the right hand side of the Carleman inequality (7). For 
the T\x\{dia l i)(dj(/))v contribution we argue analogously and include the resulting 
error contributions in the right hand side of the Carleman inequality. □ 


Remark 6. Using regularizations f m £ L p with p > n + 1 of the inhomogeneity / 
and regularizations £ W 1,p with p > n + 1 of the metric a ZJ , such that 

fm ~> f £ L n+1 (A+ r ) and -+ a lj in W 1,n+1 {A+ r ) as m —> oo, 

would have allowed us to work with C 1 solutions, w m , of the regularized problem 
(6), since on the regularized level the regularity theory of Uraltseva [Ura87] is avail¬ 
able. Passing to the limit in the regularization parameter and observing that (7) 
does not depend on the higher L p norms of the regularizations, would have pro¬ 
vided another strategy of deriving the Carleman inequality (7) and of interpreting 
the associated boundary conditions. 


Remark 7 (Antisymmetric part). It is possible to strengthen the Carleman inequal¬ 
ity (7) by including the antisymmetric part into the left hand side. As in [Rull4a], 
an additional bound can be deduced from the antisymmetric part of the Carleman 
estimate for functions which are supported in an annulus of the form Af x r . In 
Cartesian coordinates this estimate reads 
,2 


T- |||x| 1 e T<t ’w ||2 2(j 4+ r2 ) ^ (ln(r 2 /ri))“ ||e 

Thus, for functions which are supported in Af r , 
Proposition 3.1 can be strengthened to 

(23) 


2 (A+, r2 )- 

the Carleman estimate from 


e T Vr 1 (l + ln(|a:|) 2 r»t« 

+ T(lii(r 2 /ri)) _1 


L 2 (A+, r2 ) 


+ T 2 


e T ^(l + ln(|x|) 2 ) 2 Vic 


Ixl 1 e T ^w\ 


— c o 


1 C(n) (r 2 C(a ij ) ||e T 


L 2 (A+, r2 ) 


Ixl 7 1 w;| 


L 2 (A+, r ) 


|e " r|X l/|L 2 (A+, r2 )) 


p r 0 | 


for all r > 1. We will use this observation in Corollary 3.1. 


In Section 4.3, we will search for homogeneous blow-up solutions. The argument 
for the existence of these homogeneous blow-up solutions will rely on an improved 
Carleman estimate, which states that if a solution is far from homogeneous (c.f. 
(24)), it is possible to work with non-convex weight functions (c.f. condition (b)). 
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Lemma 3.2. Let t > 1, n > 2. Let p,r be two radii with 0 < p < r < 1. Assume 
that w £ with supp(w) C A^ r satisfies the following condition: There exist 

p > 0 and e £ (0,1) such that for all j £ {1,..., \ \w(r/p)\ — 2 } 

(24) \\(x-V - p)w\\ L2{A + ) >e|M| i2(A +. )• 

Moreover, suppose that it is a weak solution of the divergence form equation 

diaf 3 djW = f in A+ r , 
w > 0, d n +iw < 0, wd n +\w = 0 on A' pr , 


where f : —>■ R is an in A^ r compactly supported L 2 {A^ r ) function and 

a 13 : A+ r —> Wfyff is a tensor field such that 

(i) a i3 £ W 1 ' n+1 {A+ r ), 

(ii) the off-diagonal assumption (Al), 

(iii) the uniform ellipticity assumption (A2), 

(iv) and the following smallness condition hold: There exists S = 5(n,e,p) > 0 
such that 

t 2 sup II Va 13 II Tn+1 , ,+ ^<5 and 

p<r<r ^ 

(25) ~ 

r 2 sup | af 3 {x) — S lJ | < S. 

p<r<r 


Let (f : R —> R, </> £ C 4 , be a radial weight function satisfying the following assump¬ 
tions: 


(a) \rfi'(f) + p+ ^ < b 

(b) Tfi"(t) >~Jq, 

(c) T\$"\t)\<b 


where e := j^e ( ra +^). Set 4>{x) = </>(ln(|m|)). Then we have 

(26) ||e T 0 M _ 1 H| i2(A + r) + ll e ^ V HIL»(A+ r ) ^ 


e T *\x\ f 


L 2 (A+ r ) ' 


Here assumption (24) allows for slightly concave weight functions satisfying the 
hypotheses (a)-(c). While this permits us to strengthen the Carleman estimate by 
using (slightly) non-concave weight functions <j>, we have slight losses in the powers 
of r that appear as (polynomial) factors in the estimate (26). 


Remark 8 . • As before, estimate (26) remains valid for n = 1 with only slight 

modifications. 

• Similarly as in Remark 1, the Carleman estimate remains valid for any 
r > 0. 


Proof. We argue along the lines of the proof of the Carleman estimate of Proposition 
3.1. However, in the setting of Lemma 3.2 assumption (24) allows us to bound the 
commutator contributions, in spite of a potential (slight) concavity of the weight 
(j). To this end we have to estimate the commutator contributions 

(27) 4T((d?0)d t u,dtu) L 2 {M) + 4r 3 ((<9 2 <)>)(<9 t </>) 2 u, u) L 2 {M) - T{{dffi)u,u) L 2 (M) . 

Due to the combination of assumptions (a)-(c) and the “non-homogeneity condi¬ 
tion” (24), it is possible to absorb these contributions into the antisymmetric part 
of the operator, whenever the weight (j) becomes (slightly) concave. Indeed, after 
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the changes into conformal coordinates and after setting u = e 2 t e T ^w as well as 
t := ln(r), L 0 := \\n(r/p)\, (24) turns into 

-T<t>-?^-t fg tU _ _ n ^ 1 U _ T g t j) 


> e 


e -r <t>~—t u 


L 2 ((t+j,t+j+l)xS“ 


Thus, for j = 1,... ,L 0 - 2, 



> e 

e -T0-V* U 


L 2 x5”) 



A 2 ((t+t,t+i+l)xS") 

: j = 1 ,... ,L 0 - 2 . 


L 2 ((t-\-j,t-\-j-\- 1) xS'?) 


- T d)-lL-±t 


n — 1 




Td t <$> u 


Recalling our assumption (a), we hence obtain 
(28) 


J. 71 — 1 , 

e -r^- — tg tU 

e 

e -T^-V* U 


L 2 ((t+j,t+j+l)xS™) 2 



L 2 ((t+j,t+i+l)xS") 


i 2 ((<+it+i+i)xs;) 


Invoking the fundamental theorem of calculus to estimate 


t+j +1 

N(U) - r</>(t 2 )| < J \T(j>'(s)\ds < (p + | 

t+i 


for all fi,f 2 6 (t + j,t + j + 1), permits us to further bound the quantities in (28): 
For j = 1,. . ,,L 0 - 2 

(99) H l9tU llL 2 ((t+l,t+i+l)xS 5 ) ^ 2 e ^ + ” 1+2 ^ H U lli 2 ((i+7,t+7+l)xS") 

- 8e H U llL 2 ((t+j,t+i+l)xS”) ' 

At the boundary annuli (which correspond to j = 0 or j = Lq — 1) we do not have 
the “non-homogeneity assumption” (24) at our disposal. Yet Poincare’s inequality 
and the compact support of u imply 


(30) 


IMU 2 (m+i) ~ l|dtw|| L 2( t|t+1 ), 

\W\\L 2 (t+L 0 -l,t+L 0 +l) ~ \\9tu\\ L 2 (t+L 0 -l,t+Lo + l)- 


Using (29) as well as (30), we can proceed to control the commutator contributions 
in (27) from the Carleman inequality. In the regions in which <f> is (locally) convex, 
we argue as in the proof of Proposition 3.1. Hence, it suffices to consider the regions 
in which <f> is not convex. Let U C M be such a region. By assumption (b), we 
have 


(31) 




> (M + 


ld t u — d?d>u > t d t id t u 

L 2 (U) L 2 (U) 


The gradient contributions in (27), are immediately controlled by combining the 
estimates (31) and (29) (or respectively (31) and (30) if U contains the boundary 
annuli). For the L 2 contributions in (27) we use (29), (31) and the conditions (a)-(c) 
(or an analogous argument in which (29) is replaced by (30)). After having absorbed 
the error contributions from (27) into \ ||Ait|| i2 ( M ), the remaining 1 ||Aw|| i2 ( M ) 
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contribution then yields an L 2 lower bound (by virtue of (29), (30) and (31)). 
Moreover, an estimate similar as in (21) (without the weight) allows us to 
complement the L 2 contributions by gradient contributions. Thus, we arrive at an 
analogue of equation (16): 


i-i„ 


I LHA+ r ) 


I 


L 2 (A+ r ) — e ’ h) \\ e <l> \ x \ Aw \\mA+ r .) 


We stress that in contrast to (16), we now do not face logarithmic losses since we 
have (24) at our disposal. 

After this observation, the remainder of the proof of Lemma 3.2 essentially follows 
along the same lines as in the proof of the Carleman estimate from Proposition 3.1. 
We only consider small modifications which allow us to absorb the error contribu¬ 
tions on the right hand side of (26) into the left hand side immediately. We briefly 
comment on these modifications. Instead of estimating the differences \a — <W| in 
terms of C(a ZJ ) as in step 3 of the proof of Proposition 3.1, we use the two smallness 
conditions from (25). This becomes possible as we do not have to gain a factor |x| 7 , 
since we do not have logarithmic losses on the left hand side of (26). □ 


We will use the modified Carleman estimate from Lemma 3.2 in proving Lemma 
4.2 in Section 4.3. There we deduce a condition which ensures closeness to homo¬ 
geneous functions (along sequences of certain radii). 


3.2. Consequences of the Carleman estimate (7). We continue by deriving a 
corollary from the Carleman estimate. For us the corollary will play a fundamental 
role in the following section. In particular it will be used to deduce bounds on 
the vanishing order (c.f. Proposition 4.1, Corollary 4.1 and Lemma 4.1), the upper 
semi-continuity of the free boundary (c.f. Proposition 4.2) and doubling properties 
for solutions of (4) (c.f. Proposition 4.3). 


Corollary 3.1 (Three spheres inequality). Let be the weight function from. Propo¬ 
sition 3.1. Suppose that w : B+ —> K. is a solution of the thin obstacle problem (4) 
satisfying (AO)-(Af) and that xo £ B +, 2 . Assume that 1 < r < To < oo. Assume 
that R 0 = Ro(to, n,p, || Va* J ||j>) £ (0, |) is the radius from Remark 2 and that 
0 < rq < r 2 < r 3 < Rq and rq < |r 3 . Then there is a constant C = C(n) such that 


(32) 


r’(l + |ln(r 2 )|)- 1 e^))r 2 -Mkl| i2( ^^ (K0)) 

+ t max{ln(r 2 /ri) _1 , ln(r 3 /r 2 ) _ 1 }e T ^ ln *' T ' 2 ^V ^" 1 IMIl= ( a+ , 2r .> o)) 
< (e^'^Vr 1 M l2(< 




(^o)) 


r 2l 2r 2 V 

. e r0(ln(,3)) r -l|| u;||L2(A+ 


r 3 ,2r 3 


(so)) 


Remark 9. We stress that the constant C > 0 in Corollary 3.1 is uniform in the 

points xq € Bf. 

2 


Remark 10. If additionally to the assumptions of Corollary 3.1, the estimate (24) 
from Lemma 3.2 is satisfied, then a similar inequality as in (32) holds true for the 
weight function <f> from Lemma 3.2. More precisely, let r and cj) be in Lemma 3.2. Let 
Ro = Ro(fj,,e,T,n,p, llVa^Hi?) be such that (25) is satisfied for all p,r £ (0,I?o). 
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Then for 0 < 2 r± < r 2 < ^ < Ro, 


( 33 ) 


r0(ln(r 2 )) —1 


li2 UC 2 , 2 , 2 Uo)) 


< 


C(n^,e) (e^O) rr i |MIl 2 ( a +, 2 


Axo)) ■ 


Proof. We only prove the statements for x 3 = 0 and consider the Carleman weights, 
(f> and <f>, from Proposition 3 . 1 . Let 77 be a cut-off function supported in A 5 ri u , 3r3 / 2 
with an only radial dependence such that 


( 34 ) 


77 1 in 41.3^72,5r3/4; V 0 in IR. \ A§ ri /4 j3t - 3 /2? 

2 2 
IV77I < — in A 5ri/4t3ri/2l |V77I < — in A 5r 3/4,37-3/2- 


Then, as a consequence of the radial dependence of 77 and the support assumptions 
on it, w = wrj is supported in A 5 ri /4 3r3 / 2 and satisfies 


di(a^djw) = di(a^ dj(rjw)) 

= 2 a l;i (djrf)diW + a lj {dijr])w + ( dia lj )(djri)w in 2r3 , 
d n+1 w = r]d n+ iw = 0 on A ' ri2r3 \ A„,. 


We notice that since a 1 - 7 6 W 1,p and since Rq is chosen as in Remark 2 , the as¬ 
sumptions of the Carleman estimate ( 7 ) are satisfied. Inserting w into the Carle- 
man estimate and by using the antisymmetric part of the operator combined with 
Poincare’s inequality (c.f. Remark 7 ), i.e. respectively estimating 

tIIM’MI L=(A+ i!r2 ) ^ ( 1 n(r 2 /ri))||Au|| i2(A + r2) , 

T \\\ x \~ lv \\L*(A+ r2ir3 ) Z (ln(r 3 /r 2 ))||Au|| i2(A + ^ 

where v = e T ^w, 


leads to 


( 35 ) 


r(ln(r 2 /n)) ||e T0 |a;| (wv)\\ L2(A + ir2) 

+ T(ln(r 3 /r 2 ))" 1 ||e T 0 |x|- 1 (w; 77 )|| I/2(A + ^ 
e Tr ^\x\~ 1 (l + ln(|x|) 2 ) _ ^ (77777) 

e T( t , (l -(- ln(|a;|) 2 ) _ 2 v( UI7 /) 


- T 2 


+ T 2 


i2 (<,2. 3 ) 


< c 0 1 C(n) C(a lJ )r 2 e 


«W 2 


(?U?7) I 


+ \\e TV xla 


ij (df jV )w\\ L2{K + We^x^Hd^mw)^^ 


i2 (<,2r 3 ) 

T <t>\rlniit 


r 1 - 2r 3 


| e T ^|x| [didf 3 ) (^77)77; 11 i2(A + ^ 


where, in the notation of Proposition 3 . 1 , 7 := 1 — r -^A is the Holder differ¬ 
entiability modulus of a l °. Using the fact that 0 < r < tq is bounded, we 
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may immediately absorb the first contribution on the right hand side by choos¬ 
ing i?o = -Ro(to, 7 , n, ||Va lJ || iP ) > 0 sufficiently small as in Remark 2. Using the 
properties (34) of 77 and the monotonicity of (j>, we infer 

1 „ ^ n _j_ 

T'i 

\e T ^Vw\\ 


- n ll e IIl 2 (a; 


5t"i / 4 , 3 ri /2 


) + r ll e W \\ L 2 (A+ 


(36) 


l eT0v HL 2 (A+ 


^ ^ /2 ^5r 3 /4,3r 3 /2) 


5r 3 /4,3r 3 /2) 

_L 

5r 1 /4,3r 1 /2) 

< —e T ^ ln(ri )) ||u)|| + + l e ^(ln(r 3 )) ll^ll 

n ^ ^ ori/4,Sri/2^ 7*3 ^ V A 5r 3 /4,3r 3 /2 ■> 

- e^nUd) ||V W || i2(A+ } + e^nUa)) || Vw || i2(A+ } . 

' 5ri /4.3ri /2' ^ 5rQ/4.3rQ/2' 


br\/&,3r\/2' ' 5r 3 /4,3r 3 /2' 

In order to control the last term in (35), we use the support condition on 77 to pull 
out the exponential factor e TI ^ and to estimate the remainder via a combination of 
Holder’s inequality and Sobolev’s inequality: 

(37) 


- Il^ a "'!lrn+1^+ JMz|(0j»7)|| 2 ("+d 


\ ^ I Pi ® 1 ^ || / J rl + 1 ('A + I 

^ri,2r 3 -' v. ri,2r’ 3 -' 


i "- 1 Wi.arg) 


|||x|(9 i a y )(a j 77)u;|| L2(A + 

< C(n) ||fto y || i » +1(A +(l|(V 2 ^kh|| L2( A+ 2r3 ) + II wV7 ?Ht 2 (a+ , 2r3 ) 


+ IINVW Vr/|| L2(A + 2rg) ) • 


Due to the support assumptions on 77 and the choice of Rq, this results in estimates 
of the type (36). Pulling out the exponential factors in the remaining contributions, 
recalling the explicit form of the Carleman weight and combining (36) and (37) with 
Caccioppoli’s inequality (for which we use the boundary conditions for w and 77 on 
B[), e.g. 


11Hi 2 ! Ad 


-1 


5ri /4,3ri /2 


U 2 (A; 




results in 

3 


r ,(l + |ln(r 2 )|)- 1 e^ ln ( r2 »||7u|| i2( A+ 


r 2 > 2r 2 ' 


(38) +Tmax{ln(r 2 /r 1 )- 1 ,ln(r 3 /r 2 )- 1 }e^ 1 n(r 2 )) r -i 


+ ) 
r 2 > 2r 2' 


< Cg 1 C(n) (e^^Wrf 1 |MUa+ 1 + e ^0n(n.)) r -i | H | 


U 2 (A+, 2ri 


L 2 (A 


r 3> 2r 3 ' 


This is the desired estimate. 


□ 


3.3. Proof of Lemma 3.1. We conclude this section with the proof of Lemma 
3.1. Here we closely follow the ideas of Uraltseva [Ura87], however carefully keeping 
track of the powers of r which are involved. 

Proof of Lemma 3.1. Let ( be a smooth radial function supported in i?i( 0 )\i?i/ 8 ( 0 ), 
idendically 1 on B 1 / 2 (0)\B 1/ / 4 (0). We multiply (13) by ( 2 (2 m ~ 1 x)v and integrate 
to infer 


1 V(C(2 m - 1 >)|| i 2 <r|| 


-’IU 2 (b 2 


L \ B 2 _ m _ 2 


IIC 2 /II 


L 2 - 



THE VARIABLE COEFFICIENT THIN OBSTACLE PROBLEM 


21 


Then standard arguments show that we obtain the desired estimate (14) if we prove 
it with an additional gradient term on the right hand side of (14). 


Step 1: Energy estimates. We consider the following penalized problem 

diOgdjV* = g e in B+, 

(39) d n+1 v e = f3 e (v £ ) on B[, 

v e = 0 on dB+ . 

Here /3 e £ C°°(K) is chosen such that 

Pe <0, /?( > 0 , P e (s) = 0 for s > 0, /3 e (s) = e + - for s < —2e 2 . 

The function a\ 3 is a regularization of a lJ and g e is a regularization of 

g ■■= f - T 2 (di(f>)a lj {dj4>)v + ra lj (di<j>)djV + r{dj(j))di(a lj v) + ra 13 (dfjfyv. 

As a' e 3 is uniformly elliptic, a solution v e of (39) exists and is unique and smooth. 
Moreover, v e converges weakly in H 1 to a weak solution, u, of 

dia l3 djU = g in B+ , 

(40) dn+1 U <0, u > 0, u(d n +iu) = 0 on B[, 

u = 0 on dB+ . 

Here the boundary conditions in (40) are interpreted distributionally, i.e. similarly 
as in Step 1 of the proof of Proposition 3.1. As the bilinear form associated with 
(40) is coercive, the problem (40) has a unique weak solution. Since v satisfies (40), 
this implies that v = it, so in particular, v e — 1 v in i? 1 (l? 1 t '). 


We now proceed to the H 2 estimate for v e . Multiplying the equation with a 
compactly supported test function 77, we start with the identity 


a 13 djV^dipdx 


- J /3e{v e )r/dU n 1 - J g e gdx , 

B i B+ 


where, as above, g e is a regularization of g: 

g = f — T 2 {dicj^a 13 (dj</))v + tcl 13 {di<j))djV + r{dj<j))di{a 13 v) + ra 13 {d 2 ^(f>)v. 

We will prove 

H V 2 v£ |lL 2 (B+_ m \B+_ m _ 1 ) - Ct2 HM 2 ^HL 2 (B+_ m VB+_ m _ 2 ) 

+ c ll5 e |lL2(B 2 +_ m \B+_ m _ 2 )' 

Choosing rj = d^fj, with p £ {1,..., n} being a tangential direction, and integrat¬ 
ing by parts, results in 

J a^dj^difjdx + j (d ll a l3 )d :jl jV £ d l f)dx = - j fi' t (v t )d^v t; qd'H n ~ l 

B+ B+ B[ 
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We then set fj := ( 2 d fjV e where £ is a smooth cut-off function which is essentially 
localized in A m for some arbitrary but fixed value of m 6 N: 

J a'idj^difiV^t^dx + 2 J a^dj^d^^d^dx + 2 J ( d fl a !J ) djv e d ll: v f QdiQdx 


BT 

+ 


J (d^a l3 )d j v e di ll v e (' 2 dx = - 


j p' f (v e )(d IJ v e ) 2 ( 2 dM n 1 + J g^d^fjdx 


B? 


B[ 


B? 


< J g e d^fjdx, 

Bt 

for which we noticed that 


J P' e (v e )(d»v e ) 


2 C 2 dn n - 1 < o. 


Using the uniform cllipticity of a 1 - 7 and Young’s inequality, we thus obtain 


J l<V( 


v e \ 2 C, 2 dx < — 2 J a‘ : ‘dnjV € d ll v e (di(dx — 2 


J {d^)d. 


j v t d l _ l v e (di(dx 


B7 


J (d^a lJ )d j v e d ill v e C 2 dx + J g e d^f)dx 


B7 


B7 


<S J \dj IJ ,v t \ 2 C 2 dx + C(S) Jwm>d* 

B+ B+ 

+ 2 J \d l j,a lJ \\djV e \\d ll v e \\C\\diC\dx + J \d^a l3 \\djV e \\d itJ ,v% 2 dx 


g e d^rjdx. 


b? 


Recalling the equation for v e , this control of mixed tangential-tangential and tangential- 
normal second order derivatives can be upgraded to a full second order derivative 
control at the expense of additional error contributions: 

J \\7' 2 v e \ 2 ( 2 dx < S J \d j)l v e \ 2 C 2 dx + C(d) J |9 M u e | 2 |5iC| 2 dx 

B+ B+ B+ 

(41) +2 J \d ll a i i\\d j v*\\d^\\(\\d i (\dx+ j \d^\\dy\\d^\( 2 dx 

B+ B+ 

+ C J Q 2 \dia % i\ 2 \djV e \ 2 dx + J g e d fJ ,rjdx. 


Step 2: Error bounds. Via Holder’s inequality and interpolation, we now address 
the error contributions which involve derivatives of a lJ and which are not contained 
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in g. Using n > 2 and the smallness assumption (5), we, for instance, have: 


b: 


Id^WdyWdi^lC 2 dx < \\dy j \\ Ln+1{A + ) \\(d^v e \\ L 2 {A +)\\(d j v e \\ Hn±i) 

L n ^ (y4 m ) 

<CT||C^^|| i2(A + ) ||V(C9,^)|| i2(A + ) . 


Absorbing the second order contribution in the right hand side, allows us to deal 
with this and similar terms. We bound the terms involved in f \g e \\d^fj\dx: By 

Bt 

Holder’s inequality 


[ \g e \\duV\dx < C{6) f 

\g e \ 2 ( 2 dx + S f 

C 2 \df_ l rj\ 2 dx 

J J 

Bt Bt 

J 

Bt 


<C(S) f 

\g e \ 2 <t 2 dx + S [ 

C\dy\ 2 dx + 5 [ |(VC)^u e | 2 dx 

J 

Bt 

J 

Bt 

J 

Bt 


Absorbing the second order derivative contribution into the left hand side of (41), 
we arrive at 


IIC|VV||| l2(Bi+) < U(||(VC)9jU e || i2 ( B +) + IICdjU e || i2 ( B +) + ||n e || i2 ( B +j 

We now pass to the limit e —» 0. Using the weak lower semi-continuity of the L 2 
norm and the L 2 boundedness of g , we obtain 


(42) 


j \V 2 v\ 2 C 2 dx < C(p,n){\\{VC)djv\\ 2 L2 

b t 

+ Kd 3 v\\ 2 L * + Ml*) + J g 2 C 2 dx. 

Bt 


Hence, it remains to bound ||<?CIIi 2 ■ H ere we make use of the bound for (j) and its 
derivatives in the dyadic annuli A m to which we can always localize by choosing £ 
appropriately. In a similar spirit as before, interpolation arguments and Holder’s 
inequality permit us to control f B + g 2 C, 2 by 


J g 2 ( 2 dx<T 2 2 2m J C 2 |Vu| 

B+ B+ 


*dx + t 2 


4r>4 m 


J C 2 v 2 dx + J 


Cfdx 


As above the second order error contributions which arise on the right hand side 
can be absorbed in the left hand side of (42) 


Step 4 : Conclusion. Using the weak lower semi-continuity of the L 2 norm, we 
obtain the desired H 2 estimates for v. □ 
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4. Consequences of the Carleman Estimate: Compactness, 
Homogeneity and an Almost Optimal Regularity Result 

In this section we present the consequences of the Carleman estimate, Proposi¬ 
tion 3.1, and of Corollary 3.1: These include the independence (of the chosen subse¬ 
quence) of the rate of vanishing at free boundary points (c.f. Proposition 4.1), the 
compactness of L 2 -rescaled (sub)sequences of solutions of the thin obstacle problem 
(c.f. Proposition 4.4), the homogeneity of the blow-up solutions along particular 
subsequences (c.f. Proposition 4.5) and an almost optimal regularity result (c.f. 
Proposition 4.8). The independence of the vanishing order on the subsequence 
along which the blow-up solution is considered, combined with the homogeneity 
of particular blow-up solutions allows to classify the lowest blow-up profile and to 
obtain growth estimates (c.f. Proposition 4.7 and Corollary 4.2). 

Apart from Proposition 4.5 and Lemma 4.2, all the results of this section are 
entirely based on Corollary 3.1. Only for the proof of Proposition 4.5 do we have 
to invoke the full strength of the modified Carleman estimate from Lemma 3.2. 

In the whole section, we assume that the assumptions (A0)-(A4) are satisfied; 
in particular, all solutions of the thin obstacle problem (4) are assumed to be (L 2 ) 
normalized and the uniformly elliptic, symmetric tensor field a*- 7 is in W 1,p with 
some p € (n + l,oo]. This ensures that the assumptions of both the Carleman 
estimate (7) and of Corollary 3.1 are fulfilled throughout the section. 

4.1. The order of vanishing and growth estimates. As a direct consequence 
of the Corollary 3.1, we obtain information on the vanishing order of solutions of 
(4). Here we use the following definition: 

Definition 4.1 (Vanishing order). For any u € L 2 (B^(x o)), S > 0, we define the 
vanishing order of u at xg as 



With the aid of the vanishing order it is possible to quantify the growth of 
solutions of the thin obstacle problem. We first prove that the limsup in the 
definition of the vanishing order is in fact a limit. 


Proposition 4.1. Let w be a solution of (4) in B+ and Xg £ B[. Then 



Proof. Without loss of generality, we only show the result for Xg = 0 and derive it 
as an immediate consequence of (32) in Corollary 3.1. Indeed, if Kg < oo, for any 
e > 0 fixed, we set r = (Kg + (n — l)/2 — e)/(l + cgit/2) in (32). Note that for the 
chosen r, 


( 43 ) 


rcb(t) , n — 1 . 

- —y — (auq H” — 77 — — 6j as t —y — oo . 
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Considering the limsup ri _> 0+ on both sides of (32) and recalling our definition of 
the vanishing order, the right hand side of (32) is bounded. Thus, 

e T 4 >( inr 2 ) r -i(i _|_ | l nr2 |)- 1 || w || i2 ^ + ^ ^ < C for any 0 < r 2 <C r^. 

Hence, taking the logarithm and dividing by lnr 2 , we have 

ln IMlL2(A+ |2r2 ) > T(f>(lnr 2 ) ln(l + | lnr 2 |) In C 

In r 2 — In r 2 In r 2 In r 2 

By (43), 


lim inf 

T2—>•()+ 


HMI LHAt 2 , 2r2 ) 

In r 2 


> Kq + 


n+ 1 
2 


— e. 


Since e is arbitrary, the above inequality implies the desired estimate. 

For no = oo, we set r = M, for an arbitrarily large M, in (32) and argue 
similarly. □ 


Remark 11 (Growth bounds). • The preceding proposition immediately leads 

to the following growth estimate: For each xq £ B[ and each e > 0 there 
exists a radius r e = r t (w) > 0 such that 

r K« 0 +T-+ £ < ||tu|| L 2 ( j4 + < r K,x o + ~^~~ e for all 0 < r < r e . 

By virtue of the L 2 — L°° estimates, this also transfers to L°° norms over 
balls. 

• We emphasize that these growth bounds are asymptotic bounds as r —> 0. 
The radii r e > 0 may vary from free boundary point to free boundary point. 
Hence, these are non-uniform estimates in the free boundary points. In 
Lemma 4.1, we will indicate how to obtain related uniform upper bounds. 

• While it is possible to obtain uniform (in the radius r) upper bounds on 
the growth of solutions (c.f. Lemma 4.1), the lower bounds can only be 
improved slightly (c.f. Corollary 4.1). 

Remark 12 (Vanishing order). For points Xq with n xo > — (which guarantees 
integrability at zero), it would have been possible to consider a related notion of 
the vanishing order in which one considers the L 2 norms on (half) balls instead of 
(half) annuli: 

r ln ( 4 +(, 0 ) u2 ) 1/2 

k Xo := lim sup-^-. 

r->o+ m r 

However, at first sight this definition is less suited for our Carleman based argu¬ 
ments. Yet we remark that the two notions of vanishing order coincide. Indeed, by 
Remark 11 given any e > 0, for 0 < r < r e = r e (w) we have 


r 2K XQ +e+n +1 2 -i(2«XQ+e+n+l) < || u || 


1=0 


- H u 'llL2(B+(xo)) 


= £i 

1=0 


I L 2 (A+_ 


(j + l) ri 2-J 


- 4 .) 


< r 2K,x 0 —e+n+l 'y ' 2 - !( 2K *o -e + n + 1 ) 
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If — < k Xo < oo, then the geometric series converge, thus k xo = k Xo ; if k Xo = oo, 

then k Xq = oo. On the other side, if — < jj Io < oo, then it is not hard to see 

that there is a sequence rj —> 0 such that 


1,1 (%■., w f)' 

ln (o) 


< C < oo. 


Then by Proposition 4.1, k Xo is finite. It then necessarily equals to k X q . If R Xo = oo, 
then k Xo = oo. This follows as otherwise n Xo < oo and hence k Xo < oo, which is a 
contradiction. 


Next we show the upper semi-continuity of the map B[ E5 x <—> k, x , as a conse¬ 
quence of equation (32) in Corollary 3.1. 

Proposition 4.2 (Upper semi-continuity). Assume that w is a solution of (f). Let 
k x be the vanishing order of w at x £ B[. Then the mapping B[ 9 x i—>• k x is upper 
semi-continuous. 


Proof. We only show the upper semi-continuity at the origin. For any other point 
the argument is analogous. We show that for any e > 0 there exists S > 0 such 
that k x < Ko + e for all \x\ < 5. The statement is trivial if kq = oo, thus we assume 
kq < oo (kq > —oo follows immediately as w £ L 2 (B+)). 

Using the notation of Corollary 3.1, we choose r = (k 0 + + e)/(l + cqtt/2). 

Moreover, we consider radii ri,r 2 ,r 3 such that 

(a) V 3 = Rq/2, where Rq is the radius from Corollary 3.1, 

(b) r 2 = r 2 (e, Ko, n , Xq) > 0 is a fixed radius, which is chosen so small that 

e T&in( 2r3 ))( 2r2 )—i(j + (i n (2r 2 )) 2 )" 1 > r“ (K0+2 ^ + ^ } , 


and that for sufficiently large C = C(n,p ) (which is possible by Remark 11) 




' 9 


(44) 


HW+ I ^ 


> c!_2 _ 

^0~t~ ~t~ e 


(c) 0 < ri <C r 2 . 

Let x £ , 4 . Applying Corollary 3.1 with xq = x, the implies 


— {{n.+ l)/2+K o + ^f) 




< r -((«+l)/2-|-Ko+e) || || 

- r l " W "L\Af i 2 ri {x)) 


+ r 


IHIl 2 (A+_2, 3 (U) ' 


Noticing that A+ /a>ra (0) C A+ /i>2ra (x) for any x £ B + we switch from an 

annulus centered at x to an annulus centered at 0: 




(45) 


^ 2 (< /2 , r2 (0)) 


< —((n+l)/2+K 0 +e) | 


— ((n+l)/2+«o+e) 
3 


'L 2 (A+ 1>2ri (x)) 


Ml 


3 i 2 T 3 (x))- 


Due to the choice of r 2 in (44) and the L 2 normalization of w, the second term 
on the right hand side of (45) can be absorbed in the left hand side of (45). As a 
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consequence, we obtain 

— (( r &+l)/2+K'0 + ^f) | 


w \\lha+ 


r 2/ 2 > r 2 


(o)) 


< 


— ((n+l)/2+«o+e) | 


HI L*{A+ 




oo) 


which using (44) - we read as a lower bound for the right hand side (with r 2 
fixed): 


(46) 


C < 


— ((n+l)/2+Ko+e) 




Defining <5 = 72/4, taking the logarithm of both sides of (46), dividing by ln(ri) 
and passing to the limit r\ —> 0, yields the statement of the Proposition. □ 


The upper semi-continuity of the vanishing order immediately implies a lower 
bound for solutions around free boundary points: 

Corollary 4.1 (Growth estimates-lower bound). Let w be a solution of (f) in B+. 
Let K C B[ be a compact set. Suppose that for some R > 0, finite, n x < R for all 
x £ K. Then there exist a neighborhood, U, of K and a radius ro = ro(K, R,w) > 0 
such that 


sup \w\ > r K for all 0 < r < ro and all x £ U. 

B+(x) 


Proof. By the upper semi-continuity of the vanishing order (Proposition 4.2), there 
exists a neighborhood U C B[ of K with U compact, such that n x < R for any 
x £ U. Then by Remark 12 and by arguing similarly as in Proposition 4.2, we 
have: For any x £ U, there exists <5 = S(x, R,w) > 0 such that 

\ 1/2 


i 


B?(x) 


> r K for any x G B' s (x) and 0 < r < 6. 


By compactness of U, there exists ro = ro(K, n, w) > 0 such that 

\ 1/2 


i 


B+(x) 


> r K for any x G U and 0 < r < rg. 


Since sup B +^^ |w| > (f B +M w 2 ) 1 ' 2 , the conclusion follows. 


□ 


Remark 13. We emphasize that the lower bound in Corollary 4.1 is a non-uniform 
estimate in the solution w. In the next lemma we show a uniform upper bound. 

Lemma 4.1 (Uniform upper growth bounds). Let w be a solution of (f) in B^f. 
Given a finite constant R > 0, there exists a constant C = C(R,n,p, || Va I /|| iP ^ B +- ) ) 
such that for all xo € Tu, D B ±/2 

sup \w\ < CV min{K *o>'U| in( r )| 2 for all 0 < r < R 0 , 

B+(xo) 

where Rq = Ro(R,n,p, || VaU || iP ( B +)) is the radius in Corollary 3.1 with tq = R + n. 

Proof. We assume that 0 G P™ and will show the statement at 0. For the other 
free boundary points it follows analogously. 

Case 1: kq < R. By Proposition 4.1 for each e > 0 
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This implies that given a sequence {ej}, ej —> 0, there exist corresponding radii 


0 such that 


(47) 


-(k,o+ n ^ 1 ~ e j) 


W 


l2(a+ ) - 1 f° r all J S N. 


Let 0 < r < i?o/4- We apply Corollary 3.1 with the radii r 2 = r, r\ = Tj , with 
j 6 N such that 0 < rj -C r, and r% = Rq/2. As the weighting factor we choose 

( n— 1 \ / 1 ln(l + ln(r,) 2 )\ 1 

(48) Tj = H - --e, J |^1 - Co arctan(ln rj ) + -c 0 -^- J , 


i.e. Tj is such that e Ti ^ larj ' ) rj 1 = r 3 


— («0 + Ii 5 : i—o) 


where Cj is associated with rj as 


discussed in (47). Here Co > 0 is the constant from the Carleman inequality from 
Proposition 3.1 which will be determined later. With these assumptions, Corollary 
3.1 then yields 

e DA0nM) r -l (1 + ln 2 (r)) -i|| w || i2( ^ ) 


< 


C^C(H,p) (l + e r^(ln(Ro/2)) i? - 1 || w; | U2 (^ o/2 Ro )) . 


We note that Tj («o + 2i ^4)/(l + cot/ 2) as j —> oo. Passing to the limit j —> oo 
and using the explicit expression of (j), implies the upper bound 

| n-l \ l-c Q arctan(ln(r)) 1 1 <=0 («Q + (w -1)/2) 

r V K0+ 2 ) i+=o”/ 2 1 (l + ln 2 (r)) 2 2 r+So^T 2 ||iuj 

(49) 

< c 0 C(n,p) 

The asymptotics of the function t i-a arctan(t) yield 


LHA + 2r ) 


l + i?o (KO+2 ^ i) | 


I L 2 (A 


i? o/ 2 ’- R o 


)■ 


1 — cq arctan(ln(r)) 


1 + cof 


= 1 + 


1 


1 


1 + co7t/ 2 ln(r) 


ln(r) / 


Considering Co = (k + (n — l)/2) 1 , we then improve (49) to 

r -(« 0 +2±i)^ + ln ( r ) 2 )-i || u ,|| i2(A + 2r) 

< C(n,p,/c) ( 1 + i? 0 " (KO+ ~ ) |klli 2 (A + o/2 Jio )) • 

Recalling the parameter dependence of Rq and our assumption (AO), we obtain a 
constant C = C(R,n,p, || Va y || iP ( B +)) such that 

r ~{ K o+-r- )(i + ln(r) 2 )" 1 |MIl 2 (a+ 2 ) < C for all 0 < r < R 0 , 
or in other words 

(50) r~~^~ IMIl 2 (a+ 2 ) — ^r K ° | ln(r)| 2 for all 0 < r < R 0 . 

This then amounts to 


sup|w;| < CV K °| ln(r)| 2 , 

B+ 

by arguing similarly as in Remark 12 and using the L 2 — L°° estimate of the solution, 
i.e. sup B + |to| < C'r _( ” +1 ^ 2 ||u;|| L2 ( B +^ (c.f. [Ura87]). 

Case 2: kq > R. We define Tj := (k + ^^-4 — tj ) (1 — Co arctan(ln(r :? ))) 1 in (48) 
and argue similarly as in case 1. □ 
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4.2. Doubling and the blow-up procedure. The next proposition is a central 
result of this section. It states that at a free boundary point with finite vanishing 
order, a solution, w, of (4) satisfies an L 2 -doubling estimate. This implies sufficient 
compactness properties in order to carry out a blow-up procedure. 

Proposition 4.3. Let w : B+ —> R be a solution of (f). Then for all Xq £ r,„nBjy 2 
with k Xo < oo, there exist a radius r Xo > 0 and a constant C Xo > 0 depending on 
Xo,') ,K Xo ,n, ||Va iJ || iP ( B +) such that for all 0 < r < r Xo 


/ 


iv 2 dx < C Xo 


/ 


1 dx. 


Bf r {x 0 ) 


B+{x 0 ) 


Remark 14. We stress that, as we are using the growth estimate from Corollary 
4.1, the doubling constant C Xo > 0 and the radius r Xo > 0 are not uniform in the 
respective free boundary point. Moreover, they strongly depend on the respective 
solution w. 

We prove the doubling inequality as a consequence of Corollary 3.1. 

Proof. Without loss of generality we prove the statement of Proposition 4.3 for 
x 0 = 0. We start from Corollary 3.1 (using the second left hand side term with 
ri ~ T 2 ) and choose e > 0 such that the growth estimate from Corollary 4.1 implies 

(51) 


llw|| r2 ,o+> > e"s 1+Ko+ 2 . 


I L 2 (B+) ^ 

Moreover, we set r = + kq + 1. Then, using the notation from Corollary 3.1, 

we obtain 


'b 2 (a+) 


< C(< 


d^+Ko+lMM^)) -1 


'B 2 (A+ 2 J 


+ e(^ +K0+ m H r 3 )) r -i\\w\\ ). 

V r 3 ,2r 3 ' 


Adding 2 + K o+i)0(2e) e 1 ||w|| i 2 ( B +) to both sides of the inequality and invoking 

the monotonicity of f> yields 

(52) 


e (^i+Ko + l) 0 (ln( 2 e)) e -l || w | 


L 2 (B+) 


< C(< 


,(^i+ K0 + l)*(ln(e/2)),-l 


II ui IIl 2 (b+) 




-+K O +l)0(ln(r 3 )) 7 ,-l ll^ll 


r 3 ,2r 3 


,)■ 


Using that, by the choice of e, the growth estimate (51) holds, and potentially 

4 4 f ^ ^ j ^ j -j^ \ 

making e smaller by requiring 0 < e < Cmin{||ui|| i2( , B+ ^ , l}r 3 ' 2 , we note 

that by (51) 


= ( z V ; + K o + 1 ) 0 (ln(r' 3 )) -1 


HI; 


< _I_ e ( I1 2 i: +«o+l) 0 (ln( 2 e)) e -l ll^n 


'3 l |U ' l| £' 2 (AJ 3 ,2r 3 ) - 2C c . c II-IIL 2 (B+)- 

Thus, absorbing the second right hand side contribution of (52) into the inequality’s 
left hand side and using the explicit form of (f> once more, results in 

IMI L 2 (B+) — C II w IIl2( B +) ■ 

This is the desired doubling inequality for ro = e. □ 

The doubling property now provides sufficient compactness in order to carry out 
a blow-up procedure: 
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Proposition 4.4 (Blow-up limit). Let w be a (non-trivial) solution of (4) and let 
0 G with kq < oo. Consider the rescaling 


w a (x) 


j(ax) 


a ^ llwl 


L 2 (B+(0)) 


Then along a sequence {crjIjgN, 

w aj (x) —> wo(x) in L 2 (B+) as <jj —> 0, 

and wo is a weak solution of the thin obstacle problem with constant coefficients: 

Awo = 0 in B+, 
m > 0, -d n+1 w 0 > 0, w 0 (d n+1 w 0 ) = 0 on B[. 

Moreover, ||wo|| L 2 ( B +) = 1; in particular, it is not the trivial function. 


Proof. The proof of the L 2 convergence relies on a compactness argument using the 
doubling property as well as the control on the contributions on the boundary B[. 
We present the details: Using the gradient estimate, we obtain 

ll Vw ll a 2 (b+) < ~2 IMIl 2 (b+) ^ J w(a n+1,J dj W )dx, 

B 'lo 

~ ^2 \\ W \\l 2 (B+) > 

where the last line is a consequence of the doubling inequality and the boundary 
conditions of the thin obstacle problem. In effect, 

• ll u U|lz2(Bf) = 1’ 

• II^’ u ’o-|Il2( B +) < C. 

Hence, (along a not relabeled subsequence) we may pass to the limit a —> 0 and 
obtain w a — > wq strongly in L 2 via Rcllieh’s compactness theorem. Furthermore, 
as a u (0) = , wo weakly solves 

Awo = 0 in £?+, 
wo > 0, -dn+iwo > 0, wo{d n+ \wo) = 0 on B[. 

This proves the proposition. □ 


Remark 15. Similarly as in Proposition 4.4 it is possible to define blow-up limits 
at any free boundary point xq 6 T^,: 

w a , Xo {x) := -* w Xo {x) in L 2 (Bf) as er -A 0. 

2 II w IIl 2 (b+(x o) ) 

Similar as above, the corresponding blow-up limits are weak solutions of constant 
coefficient thin obstacle problems (however not necessarily with a®- 7 = 6 ® J ). 


4.3. Homogeneous blow-ups. In this section, we identify the minimal growth 
rates of solutions of (4) close to free boundary points. For this we show that at 
each free boundary point there exists a blow-up sequence such that the correspond¬ 
ing blow-up limit is homogeneous (c.f. Proposition 4.5). Here Lemma 4.2 plays an 
important role. Moreover, we show that the homogeneity is determined by the or¬ 
der of vanishing. Making use of the existence of homogeneous blow-up solutions as 
well as the classification of homogeneous global solutions (Proposition 4.6), we then 
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identify the lowest possible rate of vanishing (c.f. Proposition 4.7 and Corollary 4.2). 


The proof of the existence of homogeneous blow-ups is the only point in this 
section where we rely on the Carleman estimates from Section 3 directly (in the 
form of the modified Carleman estimate from Lemma 3.2), instead of only applying 
Corollary 3.1. We have: 


Lemma 4.2 (Almost homogeneity). Let w : B(~ —>■ R. be a solution of (4) for which 
the assumptions (AO) and (A2) are satisfied. Let kq > 0 be the vanishing order at 
0. Let Rq = Ro(no,n,p, ||Vab|| iP ) be the radius determined by (25) in Lemma 3.2. 
Then we have the following statements: 

(i) Let ro,rr be two radii with 0 < ro < rr < i?o and let R := ln(ri/ro). 
Assume that for some e > 0, 


(53) 


HI 




) A r 0 


Then if R = R(r±, e, n, no) is large enough, there exists r £ (ro,?r/2) such 
that 


II* • Vu; - K ow\\ L2{A + 2r) < e |M| i2(A + 2r) ■ 

(ii) There exists a sequence of radii rj -4 0 and ej —> 0 such that 
\\x- Vw - k 0 w\\ l2{a + j <C'HHl 2( a+ 2 ) f or aU 3 eN ' 

Proof. Proof of (i): We argue by contradiction. Let to := ln(ro), t\ := ln(ri) and 
let r := y/rWi correspond to t := ■ Note that t 0 + R = t\. 

Claim: If R = R(t±, n , e, kq) is large enough, then 

(54) \\ w \\ L2{Ahf) < f «o+^+T%. 

Proof of the claim: Consider the following weight function: 

(55) r4>{t) := - (k 0 + ^ t + ^{t - t 0 )(h - t). 

We aim at applying this as a (non-convex) Carleman weight in the annulus A+ ri . 
For this we have to ensure that the conditions of Lemma 3.2 are satisfied. As 
<p"" = 0, we only have to ensure conditions (a) and (b). Setting r = kq + 
condition (b) is satisfied, for R = R{no,n) chosen large enough. Condition (a) 
holds since 


T<t>'(t) = - ( « 0 + 


n — 1 


16 


-77 ) + ^(h-t)-—(t-t 0 ) 


2 R 


-K 0 


2R 
n — 1 


3e 

— 16’ ~ K ° ~ 


n — 1 


5e 

16 


Thus, the assumptions (a)-(c) of Lemma 3.2 are verified and Lemma 3.2 is applicable 
(in particular in the form of equation (33)). Consequently, we obtain 


|e’ r |x| H| l2(a + 2 _) < C 


|e T H| H 


|e T,i |d H 


r 0> 2r 0 ' 
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where C = C(n, e, kq). Inserting the specific form of the weight and the relation 
between ro,n (and to,t\) then results in 


e^f- 1 


(56) 


'L 2 (A+) ^ 


e -( Ko +^ i -T%)to ||«/|| + 

' rn 


rQ,2ro > 


_|_ e ( K o+ 2 is)* 1 


'l 2 (a +, 2ri ) ■ 


By Lemma 4.1, the RHS of (56) is bounded. Hence, 


w 


n + 1 i e \ + 


L2(A+ W ) < Ce~ T ^f = Ce^ 0+1 

— Ce^ KoH ^ 

— (^g( K oH 3— + § C*i — t) - 

= J - +16 e _ 8 tl_ 16 . 


Thus, if we take I? large enough such that 


„ 1 „ 16 In C' 

Ce 81 is < 1, i.e. R > -2 t x 


then we obtain (54). 

Now we apply Corollary 3.1 with t = (kq + (n — l)/2 — e/64)/(l + cqtz/2) and 
the three radii rj <C ro < f: 


e r0(to) ro !| l nro | 1 |k|| i2(j4 + 2r ) 


< C 


(ko + (^+1)/2—e/64). 


r 0 > 2r 0' 

H| i2(A + 2r .) -L ^-(-o+(-+ 1 )/ 2 - e / 64 )|| WM r2 ^ + 


L 2 (At„) J • 


In the limit r ? -A 0 the first term on the right hand side of the previous inequality 
vanishes due to Lemma 4.1. Thus, together with (54) and enlarging the exponent 
to absorb the logrithmic term, we obtain 


C^0+ 2 2 64) 


HI^(< 2 . 0 )^ Ca¬ 


using a summation argument as in Remark 12 further yields 


C^0+ 2 2 64 


'L 2 (Btr 0 ) ~ 


Recalling r = y/rori, implies 


I' u ’IIa 2 (b+ o ) — Cr 0 


_3 _e_ _' 

2 5 64 V^T 5 


< r. 


''O + ^T^ + A 


if R is sufficiently large or n is sufficiently small. This is a contradiction with (53). 

Proof of (ii): Note that for given e > 0, (53) is satisfied for all r < r e (w). Thus 
(ii) follows from a contradiction argument and (i). □ 


Using the result of Lemma 4.2 then yields the desired homogeneity result along 
certain subsequences: 

Proposition 4.5 (Homogeneous blow-ups). Let w : B+ -A R. be a solution of 
the variable coefficient thin obstacle problem (f). Let Xo € T^. Then there exists 
a sequence of radii {rj}jgN with rj —> 0 as j —> oo, which depends on xo, ie. 
Tj = rj(xo), such that the associated (L 2 -)blow-up sequence {uv^xoljeN converges 
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against a homogeneous solution, w Xo , of a constant coefficient equation of the type 
(2) (however not necessarily with a 13 = 8 l °). The homogeneity of w Xo corresponds 
to the degree of vanishing, k Xo , of w at Xq. 


Proof. Without loss of generality we assume that Xo is the origin. By Lemma 4.2 
(ii) there is a sequence —> 0 such that 

\\x- Ww-k 0 w\\ lha + } 
lim --- r -^~ = 0. 


3~> 00 


'LHA+, er .) 


As 11^11^2(^4+ ) < ||w|| i2 ( B + this condition, in particular, enforces 

\\x-Vw~K 0 w\\ L 2 iA + ) 

lim -- r -^~ = 0. 

J^o° \' W »L 2 (Bt j ) 

In terms of the blow-up sequence U) 2 rj the above equation can then be rewritten as 

(57) Um^ ||x • V U)2rj - K 0 W2rj \\ L *(A+ /3 J = °' 


By Proposition 4.4, up to a subsequence, W 2 rj —> wq in C^“(]R rl+1 ), where wq is a 
global solution to the thin obstacle problem (2). Hence, (57) implies that 

x ■ Vwq — Ko w o = 0 in Ay 2 v 

Therefore wq is homogeneous in AA, 2 1 . By analyticity (of solutions of the constant 
coefficient equation in the interior of R" +1 ), wq is homogeneous in R” +1 . □ 


Next we discuss the lowest possible homogeneities which appear in the previously 
described blow-up process. For this we recall (a slight modification of) Proposition 
9.9 in [PSU12], which characterizes global homogeneous solutions of (2). 

Proposition 4.6 (Proposition 9.9 in [PSU12]). Let wq be a homogeneous global 
solution of the thin obstacle problem (2) with kq € (1,2). Then Ko = 3/2 and 

w 0 {x) = Re(a; ri + ix n+ i) 3/2 

up to multiplication by a constant and a rotation in R ra . 


Remark 16. We remark that the condition kq > 1 is always satisfied if wq is a 
blow-up limit, i.e. if there exists a blow-up sequence w rj , where w solves (4) with 
0 € P„, such that w rj -A wq in C' i 1 /“(M” +1 ). Indeed, since 0 € r„, , we infer 

w rj (0) = \Vw rj (0)| = 0 by the Signorini condition and the C regularity of w rj . 
Here the identity \Vvi rj (0)| = 0 follows by noting that the function x' +->- w(x', 0) is 
C 1,0! (R(/ 2 ) regular and attains a minimum at x' = 0. The vanishing of the normal 
component of the gradient follows by approaching the free boundary point 0 € T^, 
from the interior of the open set f l Wr in combination with the C 1,a (B(~) regularity 

of w r ,. Thus, passing to the limit and using the C convergence of w rj to wq, 
yields wq(0) = |Vioo(0)| = 0. This together with the C 1,a regularity of wo implies 
that kq > 1 (more precisely kq > 1 + a). 

Proof. As the proof of Proposition 4.6 is well-known, we only give a sketch of it. 
Here we slightly deviate from the strategy presented in [PSU12], 

We recall that the key in the proof of Proposition 9.9 in [ , SU12] is to show that 
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for an arbitrary but fixed tangential direction, e £ S n D {x n +i = 0}, the associated 
directional derivative, d e wo, does not change its sign if kq £ (1,2). This then entails 
that any solution wq only depends on two variables, where one is a tangential 
direction in S n fl {x n +± = 0} and the other is the normal direction, a; n +i- In 
[PSU12] the sign condition for tangential derivatives is shown by making use of the 
Alt-Caffarelli-Friedman monotonicity formula [ACF84] . 

Indeed, this fact can also follow directly from the characterization of the second 
eigenspace for the Laplace-Beltrami operator on S n : Extending wo evenly and 
letting / denote the restriction of d e Wo onto S n , f is an eigenfunction (with A = 
(ko — 1)(k 0 +n — 2)) for Agn in S n \ A with / = 0 on An S'". Let 0 < Ai(S" \ A) < 
A 2 (S" \ A) < ..., be the spectrum of this operator. If d e wo changes its sign, then 
by the variational formulation of the second eigenvalue of Agn 

(58) A> A 2 (S"\A) > A 2 (S"). 

Since A 2 (S") is realized by the affine functions Xi with homogeneity k = 1, (58) 
implies that kq — 1 > n = 1, i.e. kq > 2, which yields a contradiction. As a 
consequence, we infer that if kq £ (1,2), then wo only depends on two variables, 
which, up to a rotation, we may assume to be x n and x n +i- Recalling the complete 
characterization of homogeneous two-dimensional eigenfunctions (c.f. [ 3 SU12]) of 
(2), we hence infer that wq{x) = Re(a; n +ix n + 1 ) 3 ^ 2 up to a multiplicative constant. 

□ 

Combining Proposition 4.5 and Proposition 4.6 yields the following classification 
of the blow-up limits at free boundary points with vanishing order (strictly) less 
than two: 

Proposition 4.7. Let w be a solution of (4) in B+. Let xo £ B[ be a free 
boundary point with k xo £ [0,2). Then k xo = 3/2 and there exists an (L 2 ) blow¬ 
up sequence w rjjX0 such that w rjjXQ —> w xo with w X0 (x) = W 3 / 2 (B~ 1 (xo)x). Here 
w 3 / 2 (x) = C n Re(x ■ v + ix n +\) 3 / 2 for some unit vector v £ S" -1 ; B satisfies 
A(x) = ( aij(x )) = B(x)B t {x) and maps {x n +± = 0} onto itself. 

Remark 17. It is possible to obtain Proposition 4.7 directly from Proposition 4.1 us¬ 
ing the Friedland-Hayman inequality [FH76] as follows (c.f. also [Andl3]): Without 
loss of generality let xo = 0. Choose e > 0 small enough such that k 0 + 2e £ (1, 2). 
By Corollary 4.1, there exists ro = r 0 (e,w) such that 

r Ko+e < w < r Ko-£ for any r £ (0,r 0 ). 

Then it is not hard to check that for every R > 1 and r^R < ro, one has 



This implies that every blow-up limit wq has a less than quadratic growth rate at 
infinity, i.e. 



sup |ui 0 | < 


< R K °+ 2e for any R > 1. 
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By the C 1,a estimates for solutions of (4), 

sup |Vwo| < R K °~ 1+2e for any R > 1 . 

Up to an affine change of coordinates, ( d e w ) ± with e £ ffi" are global subharmonic 
functions with less than linear growth at infinity. By the Friedland-Hayman in¬ 
equality, this implies that d e w has a sign. The remaining arguments are the same 
as in Proposition 9.9 of [PSU12]. 

Proposition 4.7 immediately implies the following corollary: 

Corollary 4.2. Letw : B+ ->R be a solution of (4) and assume that xo £ 

Then k Xo = § or k Xo > 2. 

4.4. Almost optimal regularity. In this section we follow a strategy described 
in the book by Petrosyan, Shahgholian and Uraltseva [PSU12] in order to obtain an 
almost optimal regularity result from our previous estimates: We begin by proving 
a in the free boundary points uniform (almost optimal) L°° growth estimate and 
then translate this into the desired (almost optimal) regularity result, c.f. Proposi¬ 
tion 4.8. 

As a corollary of Lemma 4.1 we obtain the following growth estimate in an 
i?o > 0 neighbourhood of the free boundary: 

Lemma 4.3 (Growth at the free-boundary). Let w be a solution of the variable coef¬ 
ficient thin obstacle problem (4) in B+ (0). Consider x £ B^ 2 with dist(x, Tu,) < 
Then 

|u;(x)| < Cdist(x, r ^) 3 ^ 2 ln(dist(x, P™)) 2 , 
where C = C{n , ||a lJ || whp ,p). 

Proof. This is an immediate consequence of the classification of the lowest possible 
vanishing order at free boundary points (Corollary 4.2), the L°° — L 2 estimate 
and the uniformity of Lemma 4.1 (applied with R = 3/2) in the free boundary 
points. □ 

With the growth estimate at hand we can infer an almost optimal regularity 
result which then also proves Proposition 1.1 from the introduction: 

Proposition 4.8 (Almost optimal regularity). Let w be a solution of (4) in B+. 

(i) If £ W 1,p with p > 2 (n + 1), then we have 

IVu/x 1 ) — Vu/t 2 )] < C^ 1 -x 2 | 1 / 2 ln(|a : 1 - x 2 \) 2 
for all x 1 , x 2 £ B^ 2 . 

(ii) If <W £ W 1,p with p £ (n + 1, 2 (n + 1)), then w £ C 1 , 7 (5)y 2 U B^ 2 ) with 

IVr^x 1 ) — Vwj(x 2 )| < C(y) jx 1 — x 2 | 7 for all x x ,x 2 £ B^ 2 , 
where 7 = 1 - ^. 

' p 

Apart from the dependence on 7 in the second estimate, the constants C > 0 in (i) 
and (ii) are functions of n, , ||a y ' || wljP ,p■ 
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Remark 18. • The statement (i) in particular implies that if £ W 1,p with 

p > 2 (n+l), then for any e £ ( 0 , 1 / 2 ) the function w is C 1,1 ^ 2 ~ e (B^ 2 UB ' 1 , 2 ) 
regular (with a logarithmic loss as e —> 0 ). 

• The constant C("f) > 0 in (ii) diverges as 7 \ 0 and as 7 

Up to the logarithmic loss, Proposition 4.8 provides the optimal regularity for 
solutions of the thin obstacle problem. This follows, as on the one hand the regu¬ 
larity cannot exceed the regularity predicted by interior regularity estimates with 
W 1,p coefficients. On the other hand w 3 / 2 {x) := Re(x n +ix n+ 1) 3 / 2 is a solution of 
the thin obstacle problem. 

Proof. We begin with the proof of (i) and follow the ideas in [PSU12]. For any 
x £ -Bjy 2 let d(x) = dist(x, r w ). Note that B d ^{x) n {x n+ \ = 0} must be fully 
contained in either {w(-, 0) = 0} or {u;(-,0) > 0}. Considering two points xi,X 2 £ 
B. jy 2 with lx 1 — x 2 | < where Rq > 0 is the radius from Lemma 4.3, we aim at 
showing 

(59) IVu^x 1 ) - Vu;(x 2 )| < C\x l - x 2 \i lndx 1 - x 2 |) 2 . 

For this purpose, we distinguish three cases: 

Case 1: d(x x ) > ■ Then w satisfies an elliptic equation with Holder coefficients in 

Bi(xi) with either Dirichlet or Neumann zero boundary conditions in the 
whole ball. Reflecting the metric, a 2 - 7 , and the solution, w, evenly or oddly 
with respect to {x n+ i = 0 } and using that the assumption a" +lj (x , ) 0 ) = 0 
for all j £ {l,...,n}, permits us to work with interior elliptic estimates 
(as the resulting metric remains Holder continuous). Hence, the desired 
estimate (59) is a direct consequence of this. In particular, the regularity 
of w is determined by the regularity of the coefficients a ZJ . In this case 
there is no logarithmic loss. 

Case 2: d(x 2 ) < d(x x ) < and lx 1 — x 2 \ > - . Then w solves an elliptic 

boundary value problem in the balls B d ^ x i^{x l ), i £ { 1 , 2 } (with either 
Dirichlet or Neumann data in the whole of B d t x i\(x l ) fl B[). Again we 
reflect both the solution and the metric evenly or oddly while preserving 
its Holder regularity. By Lemma 4.3 we obtain 

|w(x)| < Cd{x l )i ln(d(x 1 )) 2 in B d ^{x l ) for i £ {1,2}. 

As a consequence, interior gradient estimates for elliptic equations with 
Holder continuous coefficients imply 

|Vw;(x*)| < Cd(x i )^ ln(d( x 1 )) 2 for i £ { 1 , 2 }. 

Thus, 

iVu^x 1 ) — Vw(x 2 )| < |Vu>(x 1 )| + |Vw(x 2 )| < 2 Cd(x 1 )^ ln^x 1 )) 2 
< qx 1 -x 2 |hn(|x 1 -x 2 |) 2 . 

Case 3: d(x 2 ) < d(x x ) < and |x x — x 2 | < } . As before we combine the 

growth estimate from Lemma 4.3 with interior elliptic regularity estimates 
(after an appropriate reflection). In fact, from Lemma 4.3 we know 

sup |u>(x)| < Cd(x 1 )5 ln^x 1 )) 2 . 

Ati 1 )! 11 ) 
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Thus by the interior elliptic regularity estimates 

„1\ —1—a ii .11 


I^ Cd ( X ) 


< Cd(x 1 )" 2 ~ a \w(d(x 1 )) 2 , 


for a € (0,1). We chose a = \ (which is possible as p > 2(n + 1)) and use 
la ; 1 — a; 2 | < 1 (combined with the monotone decay of ln 2 (i) for 0 < t < 1 ) 
to infer 


iV^Or 1 ) - Vu;(a; 2 )| < \H\c^(B d(xl) (x 1 )) I ® 1 - 

2 

<C|a; 1 -a; 2 |5ln(|a; 1 -a; 2 |) 2 . 

This proves the claimed inequality for (i). 

For (ii) we argue analogously, noting that we can only apply interior regularity up 
to a C 1,7 threshold. In the argument of case 2, we notice that in By 2 we have 
la ; 1 — a ; 2 ) 1 / 2 lnda ; 1 — a ; 2 |) 2 < C(y)|a ; 1 — a; 2 | 7 . In case 3 we consider Holder estimates 
with a = 7 and notice that ^(a; 1 ) 2 _7 ln^a ; 1 )) 2 < C( 7 ). This completes the proof 
of Proposition 4.8. □ 
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